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ABSTRACT 


Measurements  are  presented  of  boundan.'  layers  with  embedded  vortices  and  with 
fihn  cooling  for  freestream  velocities  of  15,  and  11  m;s.  Measurements  of  a  boundarv' 
layer  with  embedded  vortex  and  and  no  film  cooling,  and  of  a  boundar}'  layer  with  film 
cooling  but  no  vortex  are  presented  for  freestream  velocity  of  15  nvs.  Plots  of  total 
velocity,  V,  streamwise  velocity,  V^,  secondary  flow  vectors,  total  pressure.  P^,  and 
streamwise  vorticity  are  presented  for  many  of  these   test  conditions. 

The  results  show  that  the  embedded  vortices  completely  dominate  the  flow  field 
in  boundary  layers  with  film  cooling.    This  is  indicated  from  the  plots  of  V,  V     and  P 
which  show  the  elTects  of  film  cooling  to  be  completely  decimated  in  the  vicinity  of  the 
vortex. 

In  order  to  conduct  this  study,  a  five-hole  pressure  probe  was  calibrated  for  pitch 
and  yaw.  The  probe  was  then  used  to  measure  five  pressures  associated  with  the  flow. 
From  these  pressures,  total  velocity  and  the  x,  y.  and  z  components  of  velocity  were 
determined. 

A  boundary  layer  profile  was  conducted  to  verify  the  calibration  of  the  pressure 
probe,  measurement  procedures,  and  velocity  computations.  The  results  show  expected 
boundary  layer  behavior  with  a  small  V^.  and  V^  component. 
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I.  INTRODUCTION 

The  increasing  need  for  greater  efiiciency  in  gas  turbine  engines  has  resulted  in 
higher  turbine  inlet  temperatures.  Consequently,  combustor  liners  and  turbine  blading 
are  subjected  to  greater  amounts  of  thermal  stress,  thermal  fatigue,  and  creep.  At 
present,  gas  turbines,  such  as  those  associated  with  military'  applications,  have  inlet 
temperatures  as  high  as  1800  -  2000  degrees  C  (3270  -  3632  F)  with  pressures  of  40 
atmospheres. 

Turbine  parts  may  be  protected  from  heat  loads  resulting  from  exposure  to  gas  at 
high  temperatures  by  using  a  coolant  within  turbine  passages  and  along  turbine 
surfaces.  Convection  cooling,  impingement  cooling,  transpiration  cooling,  and  film 
cooling  are  used  for  this  purpose.  Although  it  is  possible  for  the  cooling  medium  to  be 
a  substance  such  as  liquid  water  or  Freon  -  12,  most  gas  turbine  arrangements  utilize 
engine  air  bled  olT  from  the  compressor  and  rerouted  to  the  turbine  nozzles  and 
blading.  With  film  cooling,  compressor  air  is  ejected  from  surfaces  of  blades  and 
vanes.  The  film  coolant  then  protects  metal  surfaces  by  forming  a  protective  insulating 
film  between  the  blades  and  the  mainstream,  and  by  acting  as  a  heat  sink. 

The  How  through  a  turbine  cascade  is  extremely  complex.  Efforts  to  analytically 
model  the  flow  are  successful  only  when  the  fluid  is  considered  to  be  inviscid.  When 
the  effects  of  viscosity  are  included,  the  analysis  is  much  less  efiective  [Ref  1].  Because 
of  the  diOiculty  involved  in  analytical  representation  of  the  viscous  portions  of  cascade 
flows,  much  of  the  work  in  this  area  is  experimental  in  nature. 

Flow  visualization  studies  by  Herzig,  et  al.,  [Ref  2]  were  among  the  first  to  show 
the  complexity  of  flow  within  the  cascade.  Flow  visualization  studies  by  Langston.  et 
al..  [Ref  3]  and  Marchal  and  Sieverding,  [Ref  4]  also  show  the  detailed  development 
and  complexity  of  the  flow  through  the  turbine  cascade.  More  recently,  Sieverding  and 
Van  Den  Bosche,  [Ref  5]  have  used  color  smoke-visualization  to  study  the  evolution  of 
flow  in  cascades. 

Figure  1.1  from  [Ref  6]  shows  the  various  secondary'  flows  associated  with  flow 
in  turbine  cascades.  As  the  inlet  boundary  layer  approaches  the  blade's  leading  edge,  a 
horseshoe  vortex  is  formed.  The  point  at  which  this  formation  occurs  is  the  saddle 
point  which  is  clearly  shown  in  Figure  1.2  from  [Ref  4J.    One  leg  of  the  horseshoe 
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vortex  moves  initially  near  the  side  of  the  blade,  and  then  through  the  blade  passage 
where  it  is  referred  to  as  the  passage  vortex.  In  cascade  flow,  the  passage  vortex  is 
composed  of  1)  fluid  from  the  pressure  side  leg  of  the  horseshoe  vortex.  2)  the 
crossflow  from  the  endwall  boundar\'  layer,  and  3)  entrained  fluid  from  the  mainstream 
flow,  [Ref  7].  As  the  passage  vortex  continues  through  the  cascade,  it  is  forced  by  the 
pressure  gradient  to  the  opposite  side  of  the  passage  near  the  suction  side  of  the 
adjacent  blade.  This  shift  is  clearly  shown  in  Figure  1.3  from  [Ref  2].  The  second  leg 
of  the  horseshoe  vortex  follows  the  suction  side  of  the  blade  and  has  a  sense  of 
rotation  opposite  to  that  of  the  pressure  side  vortex.  This  vortex  moves  away  from  the 
corner  and  is  generally  believed  to  be  smaller  in  size  than  the  passage  vortex.  As  can 
be  seen  from  Figure  1.4  from  [Ref  8]  the  passage  vortex  makes  approximately  one 
rotation  as  it  passes  through  the  cascade.  In  Figure  1.1  the  number  of  revolutions  has 
been  exaggerated  for  clarity. 

The  objective  of  this  thesis  is  to  study  the  effects  of  embedded  vortices  on  a  fllm 
cooled  turbulent  boundan."  layer.  In  order  to  understand  the  eflects  that  this  complex 
flow  fleld  has  upon  heat  transfer,  vortex  characteristics  and  their  interaction  with 
surrounding  flow  must  be  understood.  In  this  study,  a  five-hole  pressure  probe  is  used 
to  measure  vortex  characteristics. 

Extensive  procedures  for  calibration  and  qualification  of  the  five-hole  pressure 
probe,  and  its  use  in  measuring  three-dimensional  flows  are  first  discussed.  Results  of 
six  different  test  are  then  given:  1)  a  baseline  measurement  in  a  developing  boundary- 
layer.  2)  a  boundar\'  layer  with  film  cooling  only  at  a  freestream  velocity  of  15m  s.  3)  a 
boundar}'  layer  with  a  single  embedded  vortex  and  without  film  cooling  at  a  freestream 
velocity  of  15  m,s.  4)  a  boundary  layer  with  embedded  vortex  and  film  cooling  at  a 
freestream  velocity  of  15  m  s,  5)  a  boundary  layer  with  an  embedded  vortex  and  no 
film  cooling  at  a  freestream  velocity  of  20  m,  s,  and  6)  a  boundary'  layer  with  embedded 
vortex  and  film  cooling  at  a  freestream  velocity  of  11  m,s.  Cases  3,  and  4  provide  the 
most  informative  results  of  the  study;  cases  5  and  6  are  incidental. 
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II.  EXPERIMENTAL  APPARATUS 

A.       WIND  TUNNEL 

The  wind  tunnel  pictured  in  Figure  2.1  is  now  described.    It  is  an  open-circuit 
blower  tunnel  used  to  provide  uniform  flow  at  the  nozzle  exit. 

1.  Description 

The  wind  tunnel  is  designated  the  NPS  Shear  Layer  Research  Facility  (SLRF) 
and  was  built  by  Aerolab.  It  is  designed  to  provide  uniform  flow  with  a  minimum 
amount  of  turbulence  intensity.  It  is  designed  with  numerous  pressure  taps  and  four 
38  X  20.3  cm  (15  X  8  in.)  access  ports  along  each  of  the  of  the  side  walls.  The  height  of 
the  top  wall  is  adjustable  to  permit  changes  in  the  pressure  gradient  along  the  length  of 
the  test  section.  Additionally,  the  top  wall  contains  numerous  instrument  ports  for  the 
measurement  of  various  flow  characteristics. 

The  air  speed  through  the  tunnel  can  be  adjusted  from  5  m/s  to  40  m  s.  The 
blower  exit  slips  into  the  inlet  end  of  the  wide-angle  difiuser  with  1.6  mm  of  clearance 
so  that  the  fan  is  isolated  from  the  body  of  the  wind  tunnel.  The  dilfuser  inlet  contains 
a  filter  and  nozzle.  The  test  section  is  3.048  m  (10  ft.)  long  and  0.6096  m  (2  ft.)  wide. 
The  top  is  fabricated  from  Lexan  sheet  (4.76  mm  thick),  continuously  sealed  with 
neoprene  along  the  edges.  The  tunnel's  bottom  wall  consists  of  one  1.2192  m  (4  ft.) 
long  and  three  0.6096  m  (2  ft.)  long  removeable  and  replaceable  sections.  These 
sections  are  all  0.6096  m  (2  ft.)  wide  and  are  sealed  with  "O"  rings  around  the  sides. 
Further  discussions  of  the  wind  tunnel  are  contained  in  [Ref  9]  and  [Ref  10:  p.  38]. 

2.  Qualification  and  Performance 

Extensive  qualification  test  of  the  Shear  Layer  Research  Facility  were 
conducted  by  Ligrani,  [Ref  11].  Results  show  that  the  variation  of  total  pressure  at 
the  exit  plane  of  the  nozzle  is  less  than  0.4°/o  at  26  m/s  and  34  ms.  VIean  velocity 
varies  less  than  0.7%  for  the  same  mean  freestream  speeds.  From  five-hole  pressure 
probe  measurements,  the  velocity  angle  deviation  is  nowhere  greater  than  about  0.6 
degrees  at  the  nozzle  exit  plane. 

Profile  measurements  of  the  moan  velocity  and  longitudinal  turbulence 
intensity  in  the  turbulent  boundan-'  layer  developing  at  20  m,  s  indicate  normal, 
spanwise  uniform  behavior.    For  this  qualification  test,  and  all  results  which  follow,  the 
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boundarv-  layer  was  tripped  near  the  exit  of  the  nozzle  with  a  1.5  nmi  high  strip  of 
tape.  Total  pressure  measurements  along  the  test  section  surface  at  the  nozzle  exit 
were  uniform  within  0.5%  indicating  spanwise  uniform  skin  friction. 

Freestream    turbulence    intensity   was    measure    to    be    0.00085    (8.5    one    - 
hundredths  of  one  percent  or  .085  precent)  at  20  m,  s  increasing  to  0,00095  at  30  m,  s. 

B.       INJECTION  SYSTEM 

Ordinarily,  the  injection  system  provides  film  coolant  at  temperatures  above 
ambient.  The  freestream  air  is  at  ambient  temperature;  therefore,  the  heat  transfer 
would  then  be  in  a  direction  opposite  to  that  which  occurs  in  gas  turbines.  The 
coolant  is  injected  into  the  boundary*  layer  through  a  single  row  of  injection  holes.  The 
injection  holes  are  scaled  relative  to  boundary  layer  thickness  to  be  similar  to  those 
used  in  current  turbine  blade  design. 

For  the  present  tests,  all  injected  air  was  at  ambient  temperatures. 

Injection  system  air  is  provided  by  an  71TD  Ingersoll-Rand  air  compressor.  The 
air  is  discharged  from  the  compressor  into  three  large  storage  tanks.  As  the  schematic, 
Figure  2.2.  shows  the  air  flows  from  the  storage  tank  through  an  adjustable  regulator, 
a  cut-olT  valve,  moisture  separator,  flow  regulator,  a  Fisher  and  Portor  rotometer  (full 
scale  9.345E-3  m  s,  19.8  SCFM,  model  10A3565A).  The  rotometer  which  controls  the 
volumetric  flow  rate,  discharges  the  film  coolant  through  a  dilTuser  and  into  the 
injection  heat  exchanger  and  plenum  chamber. 

The  heat  exchanger  and  plenum  chamber  shown  in  Figure  2.3  is  0.305  x  0.508  x 
0.457  m  (23  x  20  X  18  in.)  and  is  constructed  of  1.27  cm  (1/2  in.)  plexiglass.  Injection 
air  Hows  over  three  metal  plates  0.381  x  0.508  m  (15  x  20  in.).  The  two  lower  plates 
are  covered  by  silicon  rubber  heaters.  0.381  x  0.483  m  (15  x  19  in.),  rated  at  120  volts. 
The  heaters  are  controlled  through  a  type  136  Powerstat  variable  autotransformer. 
The  top  surface  of  the  chamber  contains  13  plexiglass  injection  tubes  each  being  8  cm 
(3.15  in.)  in  length  and  with  an  inner  diameter  of  0.95  cm  (3  8  in.).  This  corresponds 
to  a  length  to  diameter  ratio  of  8.42.  The  13  injection  holes  discharge  the  coolant  into 
the  boundary  layer  at  a  30  degree  angle.  There  is  a  three-diameter  spanwise  spacing 
between  the  center  of  each  hole. 

Further  discussion  of  the  qualification  and  performance  of  the  injection  system 
can  be  found  in  [Ref.  10:  p. 23]. 
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C.       FIVE-HOLE  PRESSURE  PROBE 

Multi-hole  pressure  probes  are  invaluable  in  the  investigation  and  measurement 
of  complex,  three-dimensional  flows.  In  particular,  the  five-hole  probe  is  well  suited 
for  measurement  of  three  mean  velocity  components  in  low  speed  incompressible  flows. 

The  five-hole  pressure  probe  used  to  measure  pressure  in  this  study  is 
manufactured  by  United  Sensors  and  Control  Corp.  (drawing  number 
DA-I25-24-F-22-CD).  The  probe  shown  in  Figure  2.4  is  0.6096  m  {2  ft.)  in  overall 
length  with  a  probe  diameter  of  0.318  cm  (0.125  in.).  It  is  constructed  of  corrosion- 
resistant,  non-magnetic  stainless  steel. 

The  five  pressure  holes  are  arranged  in  two  different  planes  which  intersect  at  the 
mutual  hole,  Pj.  The  probe  tip  is  prismatic  in  geometr>',  as  shown  in  Figure  2.4  The 
centrally  located  Pj  hole  is  normal  to  the  freestream.  The  pitch  plane  consists  of  Pi. 
P^,  and  P^  pressure  holes.  While  pressure  holes  P|,  P2,  P3  constitute  the  yaw  plane. 
The  distance  seperating  P2  and  P-^  is  0.1 78  cm  (0.070  in.),  P^  and  P^  are  0.155  cm 
(0.061  in.)  apart.  The  central  hole  P|  is  0.648  cm  (0.255  in.)  from  the  bottom  oi"  the 
probe  tip.   These  distances  were  measured  using  a  micrometer. 

For  calibration,  the  probe  was  positioned  in  a  manual  traversing  unit 
manufactured  by  United  Sensor  and  Control  Corp.,  Figure  2.5.  The  unit  was  modified 
by  the  user  to  include  a  compass  rose  with  a  radius  of  15.24  cm  (6  in.)  and  a  range  of 
yaw  angles  from  -40  to  +40  degrees  (0.25  degree  accuracy).  Additionally,  the  manual 
traversing  unit  positions  the  probe  vertically  in  the  mainstream  0  -  30.48  cm  (0  -  12  in.) 
range.  The  manual  traversing  unit  is  mounted  on  top  of  a  spanwise,  horizontal  sled. 
The  sled  is  62.23  cm  (24.5  in.)  wide  and  is  designed  to  set  in  place  on  top  of  the  wind 
tunnel's  side  walls.  The  spanwise  mounting  sled  is  designed  such  that  the  probe  can  be 
positioned  10.16  cm  (4  in.)  either  side  of  centerline  in  increments  of  0.635  cm  (0.25  in.). 
In  addition  to  spanwise  positioning,  the  spanwise  sled  could  be  rotated  through  a  range 
of  pitch  angles  from  -15  to  +15  degrees  with  an  accuracy  of  0.5  degrees. 

After  the  completion  of  the  probe  calibration,  an  automated  traversing 
mechanism.  Figure  2.6,  was  used  for  probe  positioning  while  measuring  the  pressures  in 
the  experimental  test  cases.  The  probe  is  fixed  into  the  automated  traversing 
mechanism  in  a  position  of  zero  yaw.  The  traversing  mechanism  has  two  degrees  of 
movement  which  allows  a  thorough  measurement  of  the  flow  field  to  be  conducted. 
Both  the  spanwise  and  vertical  traversing  blocks  are  mounted  on  a  20-thread  per  inch 
drive  screw  and  two  ground  steel,  case-hardened  steel  guide/support  shafts.    Each  drive 
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shaft  is  directly  coupled  to  a  SLO-SYN  type  MO92-FD310  stepping  motor.  The 
motors  are  comrolled  by  a  MITAS  Two-Axis  Motion  Controller,  Figure  2.7  The 
stepping  motors  and  the  controller  are  manufactured  by  Superior  Electric  Company. 
The  controller  directs  the  movement  the  probe  in  both  the  spanwise  amd  vertical 
directions.  The  MITAS  controller  comes  equipped  with  2K  bytes  of  memory  and  an 
MC68000,  16-bit  microprocessor  which  allows  the  user  to  program  the  start,  stop. 
duration,  speed,  acceleration  and  deceleration  of  the  stepping  motors. 

D.       DATA  ACQUISITION  SYSTEM 

The  data  acquisition  system,  shown  in  Figure  2.8  rapidly  acquires  the  voltages 
associated  with  each  pressure,  converts  each  vohage  to  pressure. 

1.  Transducers  and  Demodulators 

The  probe  is  connected  through  reinforced  plastic  tubing  to  five  Celesco  model 
LCVR  dilTerential  pressure  transducers.  These  transducers  have  a  designed  pressure 
range  of  0  -20  cm  (0  -  7.85  in.)  water  and  produce  a  15  to  45  mV  volt  output  signal. 
The  transducer  output  signal  is  converted  to  a  proportional  DC  signal  by  Celesco  CD 
lOD  carrier  demodulators.  Each  demodulator  has  a  maximum  frequency  response  of 
-3dB  at  500  Hz  and  a  maximum  out  put  noise  of  lOmV.  peak  to  peak.  Each 
transducer  carrier  demodulator  combination  was  calibrated  against  a  Meridian  1.27  cm 
(0.5  in.)  horizontal  manometer  with  an  accuracy  of  0.002  cm  (.005  in.)  of  water,  to  give 
typical  caUbration  of  approximately  1.0  volt  per  inch  of  water  differential  pressure. 

2.  Computers  and  Hardnare 

A  Hewlett-Packard  85  microcomputer  was  used  to  acquire  and  process  data 
for  the  caUbration  of  the  pressure  probe.  Configured  with  64K  bytes  of  memorv'  and  a 
single  magnetic  tape  cartridge  drive,  the  FIP-S5  was  used  to  collect,  store,  display,  and 
print  the  majority  of  the  data  required  during  the  course  of  probe  calibration. 

For  the  measurement  of  the  flow  field,  a  Hewlett-Packard  Series  300,  Model 
9836S  computer  was  dedicated  to  the  data  acquisition  process.  The  HP  9836S  is 
equipped  with  a  MC68000,  8  MHz  16  32-bit  processor.  Dual  5-1,4  inch  floppy  disk 
drives,  and  IM  bytes  of  memor>'.  A  HP  7470  two  pen  plotter  was  used  for  the  graphic 
representation  of  data. 

Each  transducer/ carrier  demodulator  combination  is  connected  directly  to  a 
HP  3498A  extender  which  is  controlled  by  a  HP-3497A  data  acquisitioncontrol  unit. 
The  HP-3497A  which  provides  precision  measurement  and  process  monitoring,  is 
equipped  with  analog  muhiplexing  and  a  digital  voltmeter  with  lyiV  sensitivity. 
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Six  software  programs  were  developed  for  use  during  the  thesis.  PROCAL 
was  developed  for  use  with  the  HP-85  and  was  used  during  the  probe  calibration  phase. 
PRSACQ.  VEL,  VELC,  PLOT,  VECTOR  were  developed  for  use  with  the  HP  9836S. 
PRSACQ  is  used  to  measure  the  pressures  in  the  various  How  fields.  VEL  and  VELC 
are  used  to  compute  the  velocity  components.  PLOT  and  VECTOR  are  used  for  the 
plotting  of  results.  A  thorough  discussion  of  the  programs  requires  an  understanding 
of  the  calibration  procedures  and  the  velocity  measurement  techniques  as  discussed  in 
Chapter  Three. 
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III.  CALIBRATION  AND  MEASUREMENT  PROCEDURES 

Prior  to  using  the  five-hole  pressure  probe  to  measure  flow  velocities,  it  is 
necessarv"  to  calibrate  the  probe  to  determine  the  dependence  of  yaw,  pitch,  static,  and 
total  coefficients  on  yaw  and  pitch  angles. 

A.  COORDINATE  SYSTEM 

A  right  hand  coordinate  system.  Figure  3.1  was  established  for  use  throughout 
the  course  of  study.  The  X-axis  is  parallel  to  the  streamwise  direction  and  is  positive  in 
the  downstream  direction.  The  Y-axis  is  in  the  vertical  plane  and  is  positive  from  the 
wind  tunnel's  bottom  wall.  The  Z-axis  is  in  the  vertical  plane  and  is  positive  from  the 
wind  tunnel's  bottom  wall.  The  Z-axis  is  in  the  spanwise  direction.  The  origin  of  the 
coordinate  system  is  located  on  the  centerline  line  of  the  bottom  wall. 

Yaw,  p.  is  defined  as  rotation  about  the  Y-axis  and  was  arbitrarily  defined  as 
positive  when  the  direction  oT  the  flow  resulted  in  pressure  P-^  being  greater  than  P2. 
This  condition  corresponds  to  positive  component  of  velocity  in  Z  direction.  The  pitch 
angle,  a.  is  defined  as  rotation  about  the  Z-axis  and  is  defined  as  positive  when  P^  is 
greater  than  P^.This  condition  results  in  a  positive  component  of  velocity  in  the  Y 
direction. 

B.  EXPERIMENTAL  PROCEDURES 

The  five-hole  pressure  probe  was  calibrated  using  the  method  described  by 
Treaster  and  Yocum,  [Ref  12]. 

To  perform  the  calibration,  the  probe  is  placed  in  the  manual  traversing  unit 
which  is  mounted  on  top  of  the  spanwise  horizontal  sled  as  described  in  Chapter  Two. 
The  probe  is  normal  to  the  freestream  when  P2  is  equal  to  P-j  which  gives  P  equal  to 
zero. 

The  probe  was  manually  fixed  at  a  predetermined  yaw  angle  and  then  rotated 
through  the  pitch  plane.  At  each  point,  the  data  acquisition  system  records  the  five 
pressures.    Four  pressure  coeflicients  are  then  calculated  by  the  PROCAL  program. 

Calibration  of  the  probe  was  conducted  over  a  range  of  yaw  angles  from  -20 
degrees  to  +  20  degrees  in  four-degree  increments.  The  pitch  angle  was  varied  from  -15 
degrees  to    +15  degrees  in  five-degree  increments.    This  provided  a  cone  of  angles 
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which  was  suHicient  for  the  Hows  to  be  studied  in  this  thesis.  Calibration  was 
performed  in  the  NPS  Shear  Layer  Tunnel  at  a  freestream  velocity  of  approximately  21 
m/s  (68.89  ft/sec).  The  probe  was  located  1.6  m  (5.3ft)  from  the  boundary  layer  trip. 
This  equates  to  a  Reynolds  number  of  1.97x10  based  on  downstream  distance.  Static 
pressure  was  obtained  from  a  static  pressure  tap  on  the  tunnel  side  wall  and  total 
pressure  from  a  Kiel  probe  inserted  into  the  flow  through  the  top  wall.  Static  and  total 
pressures  are  measured  once  per  set  of  calibration  points.  Atmospheric  pressure  is 
used  as  the  reference  pressure. 

C.       CALIBRATION  COEFFICIENTS 

To  use  the  probe  for  measurement  of  complex  flows,  it  is  necessary  to  determine 
the  flow  angles,  a  and  p,  and  the  local  static  and  total  pressures.  This  can  be  done  by 
determining  four  non-dimensional  pressure  coefllcients  over  a  range  of  angles  in  both 
the  yaw  and  pitch  planes.   The  four  calibration  coefllcients  are  defined  as: 

<=Pyaw=  {P2-P3)(Pl-P)  (eqn3.1) 


%itch=  (P4-P5)''(Pl-P)  (eqn3.2) 

CPtotal=(Pl-Ptotal)''(Pl-P)  (^q^3.3) 

CPstatic  =  (P-PstaticV(Pl-P)  (^^^  ^.4) 

where 

P  =  (P^  +  P^  +  P^+P^IA  (eqn3.5) 

To  be  of  value  in  measuring  a  flow  field,  the  calibration  coefficients  must  be  a 
function  of  flow  angle  only,  independent  of  velocity  and  repeatable.  Figure  3.2  shows 
that  the  calibration  is  both  independent  of  velocity  and  repeatable  in  the  yaw  plane. 
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As  previously  stated,  probe  calibration  was  conducted  by  fixing  the  yaw  angle 
and  var>-ing  the  pitch  angle.  Theoretically,  a  calibration  procedure  of  fixing  the  pitch 
angle  and  rotating  through  the  yaw  angles  should  provide  identical  results.  However, 
in  this  study  when  the  fixed  pitch,  varv*  yaw  method  was  attempted,  the  results  were  not 
identical.  The  trend  of  Cp  -^^j^  vs.  pitch  angle  showed  significant  scatter,  and  Cpj.Qj.^| 
vs.  pitch  angle  was  not  constant  for  each  yaw  angle. 

D.       RESULTS 

Figure  3.3  shows  the  variation  of  Cp^.^^^  with  the  yaw  angles  for  various  pitch 
angles.  The  response  is  nearly  linear  for  all  yaw  angles.  The  results  show  that  values 
of  Cp^.jj^^,  for  various  pitch  angles  collapse  on  top  of  each  other  for  yaw  angles  of -12 
degrees  to  +  8  degrees.  This  means  that  Cp^^^^^^.  is  independent  of  pitch  angle  in  this 
range.  From  -23  degrees  to  -15  degrees  and  from  +8  degrees  to  +  15  degrees,  there 
are  slight  variations  in  the  results  indicating  that  Cp^.^^^  has  a  slight  dependency  on  the 
pitch  angle.  The  fact  that  yaw  is  only  slightly  dependent  on  the  pitch  angle  allows  data 
to  be  more  easily  processed  in  determination  of  fiow  velocity. 

The  variation  of  Cp_^^j^  vs.  pitch  angle,  Figures  3.4  and  3.5,  show  that  this 
coefficient  is  dependent  on  both  yaw  and  pitch  angles.  The  trend  of  Cpp-^-^u  vs.  pitch 
angle  is  generally  linear  but  there  are  variations  from  that  linearity  for  each  yaw  angle, 
and  unlike  the  yaw  plane,  these  variations  are  not  restricted  to  any  particular  region. 

Figures  3.4  and  3.5  show  that  the  range  of  values  for  Cp  -^^j^  is  small  over  the 
range  of  a  shown  compared  to  the  variation  of  Cp.^^.^^^^,  with  p.  [Ref  12:  pp.  27-28] 
attribtues  this  to  the  types  of  surfaces  on  which  the  holes  in  the  pitch  and  yaw  plane 
are  connected.  Large  yaw  angles  result  in  one  hole  being  nearly  aligned  with  the  fiow. 
This  hole  senses  a  pressure  nearly  equal  to  the  total  pressure  of  the  freestream.  The 
other  hole  is  then  blocked  from  the  freestream  and,  consequently,  reads  a  pressure 
much  less  than  the  freestream  static  pressure.  The  holes  in  the  pitch  plane  have  a 
different  response  to  variation  in  the  pitch  angle.  When  pitched,  one  hole  reads  a 
pressure  which  is  near  the  total  pressure  of  the  freestream,  but  the  second  hole  senses  a 
pressure  which  is  greater  than  the  freestream  static  pressure.  Thus,  Cp^^^.^^^^.  is  a  much 
larger  number  than  Cpp-^^j^.  The  small  range  of  Cp  -j-^j^  increases  the  scatter  and 
uncertainty  of  pitch  angle  measurements. 

The  plot  of  Cpj^^^^i^.  vs.  pitch  angle,  Figure  3.6,  shows  that  Cp^^^^^-^  has  weak 
dependency  on  yaw  and  pitch  angles. 
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^Ptotal  ^'^^  P^^^^  ^ng^^'  Figure  3.7,  indicates  that  for  any  given  yaw  angle, 
CPtotal  ''  constant  throughout  the  range  of  pitch  angles.  Here,  Cp^Q^,^j  also  shows  a 
weak  dependency  on  yaw  and  pitch  angles. 

E.       INTERPOLATION  AND  APPLICATION 

After  the  probe  has  been  calibrated  and  the  operating  characteristics  of  the  probe 
are  known,  it  is  possible  to  determine  the  pitch  angle,  yaw  angle,  local  static  and  total 
pressures  for  any  How  field, 

The  probe  is  positioned  normal  to  the  freestream,  and  at  any  location  in  the  flow 
field,  the  five  pressures  can  be  measured.  These  five  pressures  are  now  used  to 
calculate  the  experimemal  or  local  coefficients  of  yaw,  pitch,  static  pressure,  and  total 
pressure. 

A  fifth  order  polynomial  was  fitted  to  the  average  values  of  Cp  .  .  The 
resulting  polynomial  computed  using  a  FORTRAN  program  based  on  the  least-squares 
method  is: 

Papp=--158-7.36(Mi)  +  0.135(Mi)2  +  0.304(Mj)3  +  0.009(Mj)4-0.03I(M^)5  (eqn  3.5) 

^^^^'  Papp  ^s  ^he  approximate  yaw  angle,  and  Mj  is  the  local  coefficient  of  yaw. 

Knowing  the  approximate  yaw  angle  and  the  local  Cppj^^^^  a  computerized 
interpolation  is  performed  to  determine  the  pitch  angle.  Since  the  value  of  the  pitch 
angle  is  dependent  on  the  local  Cp^^^^j^  and  on  the  yaw  angle,  it  is  necessar>'  to 
perform  a  double  interpolation.  Referring  to  Figure  3.8,  the  two  yaw  angles,  Pj  and 
P2^  which  bracket  the  approximate  yaw  angle  are  first  determined.  Next,  the  local 
^Ppitch,  -^^2a  3^^  ^hb,  ^^  located  between  known  values  of  Cp  -^^j^  from  the 
calibration  data  for  each  yaw  angle.  In  Figure  3.8,  these  values  are  designated  C^ 
^12,  ^21,  ^nd  C22  which  correspond  to  pitch  angles  from  the  calibration  data  of  aj/ 
«12,  a2i,  and  022-  The  bracketing  pitch  angle  a^^^  is  determined  by  the  following 
interpolation 

(^11  ~  '^2a)'(CH  -Ci2)  =  (ajj  -aj^);(ajj  -Qj^)  (eqn  3.6) 

which  rearranged  gives 

«,f  =  «,,  -(a,,  -«,2K(Cii  -  M2,);((Cjj  -q^))  (eqn  3.7) 
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The  other  bracketing  pitch  angle  ct2f  is  found  in  a  similar  manner  which  gives 

The  pitch  angle  for  the  flow  a^is  found  through  a  second  interpolation  which  results  in 
the  relationship 

«f  =  «lf-  («lf-  «2fX(Pl  -  Papp)''(Pl  -  P2))  (^^^  3-9) 

Because  Cp^.^^^  has  a  slight  dependency  on  pitch,  it  is  necessary  to  refine  the 
approximate  yaw  angle  once  the  local  pitch  angle  is  known.  This  is  done  with  a 
computerized  interpolation  routine  which  is  ver\'  similar  to  that  used  for  the  pitch 
angle.  As  shown  in  Figure  3.9,  the  pitch  angles  from  the  cahbration  data,  a^  and  02 
which  bracket  the  local  pitch  angle,  a^  are  first  determined.  The  local  Cp^,^^^,  VIj^  is 
then  located  between  the  values  of  Cp^.^^^.  from  the  calibration  data  for  aj  This  results 
in  C-^j  and  C-»j  which  correspond  to  yaw  angles  Pjj  and  Pp.  Then,  by  linear 
interpolation  Pjf  is  found 

(C3i-M^j,)/(C3j-C2i)  =  (((Pii-Pif)'(PH-Pi2))  (eqnS.lO) 


or 


Plf=Pll-(Pll-Pl2X(^3l-^lb)'(C3i-C2i))  (eqn3.1i) 
The  same  interpolation  is  done  for  aj  which  gives 

hr  P2I  -  (P2I  -  P22K(<^32  -  ^la)  ^<=32  "  ^22))  ^^^^  ^-  ^2) 

p2r  and  P2r  are  the  yaw  angles  which  bracket  the  yaw  angle  of  the  How.  The 
flows  yaw  angle,  Pr  can  be  found  by  a  second  interpolation  which  results  in  the 
relationship 

Pf=P2f-(P2f-PifX(«2-«f)'^«2-«p)  (eqn  3.13) 
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Since  local  total  and  static  pressures  are  not  measured  at  each  probe  location,  it 
is  necessar\'  to  perform  another  double  interpolation  to  compute  their  values.  From 
the  calibration  data,  the  values  of  Cpj-Qj-^^^  for  each  of  the  bracketing  yaw  angles  and 
pitch  angles  are  known.  In  Figure  3.10,  these  points  are  designated  Cjp  Cj2.  C-)i, 
and  C-)2-  An  interpolation  is  performed  to  determine  the  Cp^^^^j  corresponding  to  the 
local  pitch  angle.  These  points  are  designated,  C||~and  €2^  The  local  Cp^Qj-^^^,  M3,  is 
found  using  the  relationship 

(C^j--  M3);(Cif-  C,p  =  (P^  -  p^V{Pj  -  Pp  (eqn  3.14) 

which  when  rearranged  as  shown  below  gives  the  local  Cp^^j.^^  for  the  flow 

^3  =  ^lf"(^lf"^2fX(Pi-Pf)(Pi-P2^)  (eqn  3.15) 

The  local  Cp^^^^^^    M^  is  found  in  a  similar  manner.    Referring  to  Figure  3.11,  the 
relationship  for  deterniining  the  local  Cp  ^^^[[q  is  then  given  by 

M^=C3j.-(C3j.-C^j-)((pj-pj.)  (Pj-P,))  (eqn  3.16) 


F.       VELOCITY  COMPONENT  DETERMINATION 

Once  the  values  of  the  local  Cp^^^^^j  and  Cp^^^.-  have  been  determined,  the  local 
total  and  static  pressure  can  be  calculated  and,  subsequently,  the  total  velocity  at  the 
probe  tip  can  be  determined. 

The  defining  relationships  for  Cp^.^^^^  and  Cp^^^^-^  can  be  rearranged  to 
determine  ?^^^^^  and  P^^^^^^: 

Ptotal=Pi-(CPtotalXPl-P)  (eqn  3.17) 


Pstatic=P-(CPstaticKPl-P)  (eqn  3.18) 

By  using  Bernoulli's  equation,  the  magnitude  of  the  local  total  velocity  is: 

->■) 


^'=v^-(^total-Pstatic)P  (eqn3.19) 

The  three  components  of  velocity  can  now  be  determined  using  the  total  velocity 
vector  and  the  local  pitch  and  yaw  angles.  Referring  to  Figure  3.1,  these  velocity 
components  are  given  by, 

V^  =  Vcosacosp  (eqn  3.20) 


Vy  =  Vsina  (eqn  3.21) 


V    =  VcosasinP  (eqn  3.22) 


G.       SPATIAL  RESOLUTION  CORRELATION 

The  calculation  of  V.^,  may  be  influenced  by  the  local  total  velocity  gradient. 
Corrections  for  this  effect  may  be  made  using  the  following  relationship 

V"^yo^(^^''^^Xl>')  (eqn  3.23) 

Here  V^.^  is  the  uncorrected  value  ofV\,  and  V  .  is  the  value  of  V,,  corrected  for 
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spatial  resolution.  The  value  of  ly  used  was  slightly  greater  than  the  distance  between 
P^  and  P^  (0.155  cm):  a  value  of  0.200  cm  gave  constant  V.^  through  the  two- 
dimensional  boundan.'  layer. 

H.       SOFTWARE 

Six  programs  were  developed  for  use  during  this  study.  They  are  PROCAL, 
PRSACQ,  VEL  VELC,  PLOT,  and  VECTOR.  Each  program  is  written  in  BASIC. 
PROCAL  was  used  with  the  HP-85,  all  others  were  written  for  the  HP-9836S.  All  of 
the  programs  are  listed  in  Appendix  C. 

PROCAL  is  a  BASIC  program  used  for  the  calibration  of  the  pressure  probe. 
The  program  begins  by  computing  the  correction  factor  for  random  noise  associated 


with  each  transducer.  The  user  is  then  prompted  to  manually  calibrate  each  transducer 
against  a  horizontal  monometer.  Static  pressures  are  input  utilizing  a  static  pressure 
tap  on  the  side  wall  and  one  of  the  transducers.  Total  pressure  of  the  freestream  is 
input  using  a  Keil  probe  inserted  into  freestream  and  a  transducer.  The  user  is  next 
prompted  to  input  the  ambient  pressure  in  inches  of  mercur>'.  After  the  user  has 
positioned  the  probe  at  the  desired  angles  of  yaw  and  pitch,  those  angles  are  input  into 
the  program. 

The  computer  then  acquires  the  five  voltages  from  the  data  acquisition  system, 
converts  each  voltage  to  a  pressure  in  inches  of  water,  and  then  calculates  the  yaw. 
pitch,  total,  and  static  coefficients  of  pressure.  Finally,  the  yaw  angle,  pitch  angle,  and 
the  four  coefficients  are  stored  on  a  separate  file  and  printed  out  utilzing  the  HP-85's 
internal  printer. 

PRSACQ  was  used  to  acquire  the  pressures  during  the  experiment.  PRSACQ 
begins  by  prompting  the  user  for  the  number  spanwise  and  vertical  data  points  and  the 
resolution.  A  matrix  of  data  point  location  is  then  computed.  Next,  the  transducers 
are  corrected  for  noise  and  calibrated  against  the  manometer.  Freestream  static  and 
total  pressures  are  measured  and  ambient  conditions  are  input.  The  program  enters  a 
loop  which  samples  each  pressure  ten  times  per  probe  location.  The  local  Cp.^,^^^,  and 
^Ppitch  ^^^  computed.  Probe  position,  Cp^.^^^^,  and  Cpp-^^i^  P|  and  P^otal  ^^^  stored 
in  a  matrix.  At  the  end  of  the  data  collection  run.  these  values  are  read  into  a  data  file 
on  a  floppy  disk. 

VEL  is  used  to  process  the  raw  data  acquired  by  PRSACQ.  The  data  is  first  read 
into  computer  memor>'.  The  program  computes  an  approximate  Cp.^.^^^  using  a 
polynomial  fit.  Next,  double  interpolation  subroutine  is  used  to  comput  local  Cp_j^  j^. 
A  second  interpolation  subroutine  is  used  to  refine  the  value  of  Cp^.j,^^..  The  values  of 
the  Cp  -^^j^  and  Cp.^.^^  are  used  in  a  third  interpolation  subroutine  to  compute  the 
local  Cpg^jjj-^.  and  Cp^Q^-^i.  The  total  pressure  is  found  from  the  definition  of  Cpj.Qj.^j. 
The  velocity  at  the  probe  tip  is  computed,  and  the  x,  y,  and  z  components  of  velocity 
are  determined.  Probe  position  (y  and  z  coordinates),  total  pressure,  total  velocity,  and 
V     V     and  V    are  stored  in  a  matrix  and  then  read  into  a  data  file. 

VELC  corrects  the  V^.  component  for  spatial  resolution.  PLOT  is  used  to 
generate  graphs  of  streamwise  velocity,  total  velocity,  and  total  pressure.  VECTOR  is 
used  to  plot  the  secondarv*  flow  vectors. 
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IV.  EXPERIMENTAL  RESULTS 

The  study  was  conducted  in  three  parts.  The  first  was  the  measurement  of 
baseline  data  consisting  of  boundan.'  layer  profiles  in  a  two-dimensional  mean  How 
field.  The  second  was  the  measurement  of  the  boundary-  layer  characteristics  at  a 
freestream  velocity  of  15  m,  s  with  1)  film  coohng  only,  2)  with  embedded  vortex  only, 
and  3)  with  both  film  cooling  and  embedded  vortex.  The  third  part  was  the 
measurement  of  boundar\'  layer  with  vortex  at  20  m/s,  and  measurement  of  boundary- 
layer  with  vortex  and  film  cooling  at  a  freestream  velocity  of  11  m/s. 

A.  BASELINE  RESULTS 

The  baseline  boundar\'  layer  profiles  were  conducted  at  a  freestream  velocity  of 
22  m's.  For  these  tests,  the  top  wall  of  the  wind  tunnel  was  adjusted  at  20  m,  s 
freestream  velocity  so  that  a  zero  pressure  gradient  existed  within  0.15  mm  water  along 
the  length  of  the  test  section.  Profiles  were  taken  at  three  spanwise  locations  z=  +2.54 
cm.  z  =  0  cm,  and  z  =  -2.54  cm.  Figure  4.1  shows  measurements  of  the  slreamwise 
velocity,  V^^.  Figure  4.2  shows  measurements  of  V  .  Figures  4.3  and  4.4  show 
mesurements  of  the  normal  velocity,  V.^..  The  results  of  the  baseline  measurements 
indicate  behavior  expected  of  a  2-D  turbulent  boundar\'  layer,  since  the  figures  show 
mean  profiles  to  be  spanwise  uniform  for  all  three  velocity  components.  Figures  4.2 
and  4.3  show  that  the  V  and  the  corrected  V.^^,  components  are  small  and  nearly  zero 
as  would  be  expected.    Figure  4.4  is  the  plot  of  V.^,,  uncorrected  for  spatial  resolution. 

B.  15M/S  FREESTREAM  VELOCITY  RESULTS 

The  investigation  of  the  flow  field  was  conducted  by  using  the  five-hole  pressure 
probe  to  measure  pressures  at  800  points  in  a  spanwise  plane.  Data  was  taken  at  20 
different  vertical  locations,  each  having  40  spanwise  locations.  All  15  m  s 
measurements  were  taken  at  a  location  of  1.49  meters  from  the  boundary'  layer  trip  or 
0.39  meters  from  the  heat  transfer  plate  leading  edge.  The  film  cooling  cases  were 
conducted  with  injection  air  at  75%  of  full  scale  on  the  rotometer  which  corresponds 
to  a  blowing  ratio  (ratio  of  coolant  to  mass  fluxes)  of  0.50. 
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1.  Boundan'  layer  uith  film  cooling 

The  results  of  boundary  layer  with  film  cooling  only  and  no  vortex  are  shown 
in  Figures  4.5  to  4.8.  Figure  4.5  shows  contours  of  the  streamwise  velocity,  Figure  4.6 
shows  total  velocity  contours,  Figure  4.7  is  the  plot  for  total  pressure,  and  Figure  4.8 
shows  the  secondary  flow  vectors.  Away  from  the  wall,  outside  the  boundary"  layer,  the 
first  three  figures  show  spatially  uniform  behavior  for  V^,  V,  and  P^.  Secondar>-  flow 
vectors  are  ver>'  small  ever\'where  in  Figure  4.8.  Near  the  wall  deficits  of  V^,  V,  and 
P  correspond  to  locations  of  the  film  cooling  jets  which  are  located  about  every  3.0 
cm  from  tunnel  centerline. 

2.  Boundary  layer  >vith  vortex 

The  vortex  is  generated  by  using  a  half-delta  wing  which  is  3.0  cm  high  with 
7.5  cm  chordand  an  angle  of  18°.  It  is  identical  to  vortex  generator  #  2  described  by 
Joseph,  [Ref  10:  p. 76].  The  vortex  generator  was  located  at  a  spanwise  location  of 
z  =  4.79  cm  (note  that  the  direction  of  +z  in  [Ref  10:  p.75]  is  reversed).  Figure  4.9 
shows  the  streamwise  velocity  results,  Figure  4.10  is  a  plot  of  the  total  velocity 
contours.  Figure  4.11  is  a  contour  plot  of  total  pressure,  and  Figure  4.12  shows  the 
secondary'  flow  vectors.  The  contour  plots  for  V,|^,V,  and  P^  show  significant  deficits 
caused  by  the  generator  wake  which  is  rolled  up  with  the  vortex.  The  center  of  the 
wake  is  located  at  y=3.3  cm  and  z  =  -3.05  cm.  Figure  4.12  shows  that  the  vortex 
center  is  located  near  the  same  location,  with  overall  characteristics  similar  to  a 
Rankine  vortex.  Figure  4.13  shows  the  streamwise  vorticity  contours,  where  the 
vorticity  is  calculated  using 

(Of  =  {dN^^dyy  (dVyldz)  (eqn  4.1) 

The  vorticity  is  largest  near  the  vortex  center  as  expected.    The  circulation  of  the 
vortex  is  estimated  to  be  0.2708  m  /s  using  the  equation 

r  =  jco^dA,  (eqn  4.2) 

3.  Boundary  layer  with  vortex  and  film  cooling 

Figure  4.14  shows  the  streamwise  velocity  contours  for  this  case.  Figure  4.15  is 
the  total  velocity  contour  plots,  Figure  4.16  shows  the  results  for  total  pressure. 
Figure  4.17  shows  the  secondary  flow  vectors,  and  Figure  4.18  is  a  contour  plot  of 
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vorticity.  Figure  4.15.  4-16,  and  4-17  show  that  the  deficits  for  V^,  V,  and  P^  from  the 
film  cooling  are  no  longer  present  near  the  vortex.  This  result  shows  that  the  eflects  of 
film  cooling  are  decimated  by  the  vortex.  The  efTect  of  the  vortex  is  particularly 
pronounced  near  its  downwash  side.  As  for  the  previous  case,  Figure  4.18  shows  that 
vorticity  is  again  highest  near  the  vortex  center.  Figure  4.19  from  [Ref  10:  p.99]  shows 
that  the  results  of  this  case  are  consistent  with  previous  work  by  Joseph  which  shows 
high  heat  transfer  rates  on  the  downwash  side  of  the  vortex  and  low  heat  transfer  rates 
on  the  upwash  side  of  the  vortex.  The  circulation  for  this  case  is  estimated  to  be 
0.2708  m^  s. 

C.       20  M/S  AND  10  M/S  FREESTREAM  VELOCITY  RESULTS 

Figure  4.20  is  a  plot  of  streamwise  velocity  for  the  case  of  an  embedded  vortex 
with  no  film  cooling  with  a  freestream  velocity  of  20  m;s.  Figure  4.21  is  a  plot  of  the 
secondar>"  flow  vectors  for  the  same  case.  Figure  4.22  and  Figure  4.23  are  for  a 
boundarv'  layer  with  embedded  vortex  and  film  cooling  at  11  m/s  freestream  velocity. 
The  blowing  ratio  for  the  20  m  s  case  is  0.38,  and  the  blowing  ratio  at  11  ms  was  0.68. 
The  results  for  both  cases  show  trends  which  are  similar  to  those  discussed  above.  The 
results  at  a  freestream  velocity  of  11  m  s  are  believed  to  be  less  reliable  because  of 
disturbances  at  the  inlet  of  the  tunnel  which  propagated  to  the  test  section  during  the 
time  period  the  data  was  acquired. 
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V.  SUMMARY  AND  CONCLUSIONS 

VIeasurements  are  presented  of  boundary  layers  with  embedded  vortices  and  with 
film  cooling  for  freestream  velocities  of  15,  and  11  m/s.  VIeasurements  of  a  boundary 
layer  with  embedded  vortex  and  and  no  film  cooling,  and  of  a  boundar>'  layer  with  film 
cooling  but  no  vortex  are  presented  for  freestream  velocity  of  15  m,  s.  Plots  of  total 
velocity,  V,  streamwise  velocity,  V,^,  secondar>'  flow  vectors,  total  pressure,  P^,  and 
streamwise  vorticity  are  presented  for  many  of  these   test  conditions. 

The  results  show  that  the  embedded  vortices  completely  dominate  the  fiow  field 
in  boundar\'  layers  with  film  cooling.  This  is  indicated  from  the  plots  of  V,  V^,  and  P^ 
which  show  the  efiects  of  film  cooling  to  be  completely  decimated  in  the  vicinity  of  the 
vortex.  This  result  is  consistent  with  the  heat  tranfer  results  of  Joseph,  [Ref  10:  p. 54], 
which  shows  a  localized  hot  spot  at  the  wall  near  the  same  location.  Future  film 
cooling  injection-hole  arrangements  in  turbine  blades  must  be  designed  to  compensate 
for  hot  spots  due  to  the  vortices. 

In  order  to  conduct  this  study,  a  five-hole  pressure  probe  was  calibrated  for  pitch 
and  yaw.  The  probe  was  then  used  to  measure  five  pressures  associated  with  the  fiow. 
From  these  pressures,  total  velocity  and  the  x,  y,  and  z  components  of  velocity  were 
determined. 

A  boundary'  layer  profile  was  conducted  to  verify  the  calibration  of  the  pressure 
probe,  measurement  procedures,  and  velocity  computations.  The  results  show  expected 
boundary'  layer  behavior  with  a  small  V.^  and  V^  component. 

It  is  recommended  that  flow  visulization  study  of  the  interaction  betw^een  the 
vortex  and  film  cooling  be  conducted  to  enhance  the  understanding  of  this  complex 
phenomena. 
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APPENDIX  A 
FIGURES 


Fisure  1.1     Endwall  secondarv  flows. 
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V2    =160m/5 


1  .      Saddle   point 

2.      Separation  line  of  inlet  end  wall  Ijoundary  layt-r 

3a.  End  wall  separation  line  of  suction  side    leu   of 

leadmu  ed^e  vortex 
3b.   End  wall  separation  line  of  "pressure   side   U-'_: 

of  leadinu   edue   vortex"  and    or   passaue   vortex 
-4.      Start   of  suction   side  end   wall  corner   vorti-x 
5.       Low  enerey   reiiit)n 


Figure  1.2     Smoke  and  oil  visulization  for  rotor  blade. 
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Fieure  1.3     Secondarv  flow  deflection. 
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Figure  1.4     Rotation  of 


passage  and  corner  vortices  in  turbine  cascade. 


Figure  2.1     P!io:v^uraphs  of  wind  tunnel. 
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Figure  2.2     Schematic  of  injection  air  flow. 


Figure  2.3     Schematic  of  injection  plenum. 
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Figure  2.4     Photographs  of  pressure  probe. 
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Figure  2.5     Manual  traversing  device. 
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Fi2ure  2.6    Automated  traversine  mechanism. 
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Figure  2.7    Photographs  of  Two-Axis  Motion  Controller. 
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Figure  2.8     Photographs  of  data  acquistion  system. 
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Figure  3.1     Coordinate  system  for  flow  measurement. 
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Figure  3.2     Repeatabilitv  of  calibration  results  at  dilTerent 
freestream  velocities   Ibr  Cp^.„„,  vs.  yaw  angles. 
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Figure  3.3     Probe  calibration.  C 
freestrcani  velocitv 
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Figure  3.4     Probe  calibration.  Cp  •,  ^  vs.  pitch  anele  for  -20o  to  -4o  vaw, 

freestream  vtlcrciitv  of  21  m,  sT 
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Figure  3.5     Probe  calibration.  Cp  -^  u  vs.  pitch  anglefor  O"  to  +20°  vaw, 

freestream  Vfeiaditv  of  21  m/sT 
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Figure  3.6     Probe  calibration,  Cp^^^^-^  vs.  pitch  angles,  freestream  velocity  of  21  m/s. 
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Figure  3.7     Probe  calibration,  Cp.  .  i  vs.  pitch  angles,  freestream  velocity  of  21  m/s. 
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Figure  3.8     Interpolation  for  pitch  angles. 
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Figure  3.9     Interpolation  for  yaw  angles. 
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Figure  3.10     Interpolation  for  local  Cp^^^^j. 
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Figure  3.11     Interpolation  for  local  Cp  f. 
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Figure  4.1     Baseline  boundary  layer  results  for  streamwise  velocity,  V, 
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Figure  4.2     Baseline  boundarv"  layer  results  for  spanwise  velocity.V. 
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Figure  4.3     Baseline  boundary'  layer  results  for  corrected  V 
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Figure  4.4    Baseline  boundary  layer  results  for  uncorrected  V 
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Figure  4.5    Streamwise  velocity  for  boundarv  layer  without  embedded  vortex, 
75%  film  coolmg,  freestreani  velocity  15  m/s. 
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Fieure  4.6    Total  velocity  for  boundary'  layer  without  embedded  vortex 
75%  film  co'^oling,  freestream  velocity  15  m/s. 
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Figure  4.7    Total  pressure  for  boundar\'  layer  without  embedded  vortex, 
75%  film  cooling,  freestream  velocity  15  m/s. 
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Figure  4.9     Streamwise  velocity  for  boundary  layer  with  embedded  vortex, 
without  film  cooling,  freestream  velocity  15  m;s. 
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Figure  4.10    Total  velocity  for  boundary'  layer  with  embedded  vortex,  without 
film  cooling,  freestream  velocity  15  m/s. 
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Figure  4.12     Secondan*  flow  vectors  for  boundary'  layer  with  embedded  vortex, 
without  film  cooling,  treestream  velocity  15  m;s. 
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Figure  4.16    Total  pressure  for  boundar\'  layer  with  embedded  vortex 
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Figure  4.17    Secondan-'  flow  vectors  for  boundary  layer  with  embedded  vortex 
and   film  cooling,  freestream  velocity  15  m;s. 
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Figure  4.18     Vorticity  contours  for  boundarv  laver  with  embedded  vortex 
and  film  cooling,  freestream  velocity  15  m,'s. 
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Figure  .4. 2 1     Secondar\'  flow  vectors  for  embedded  vortex 
without  filmcooling,  freestream  velocity  20  m;s. 
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Figure  4.22     Streamwise  velocity  for  embedded  vortex 
with  film  cooling,  frecstrcam  velocity  1 1  m/s. 
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Figure  4.23     Secondarv  flow  vectors  for  embedded  vortex 
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APPENDIX  C 
SOFTWARE 


I'j    R E 1 -^  F R G G K H M  PR 0 C h L 

2y  Rr!;'-:  THIS  PROGRRM  COMPUTE? 

7-3  REM  ""HE  COEFFICIENTS  OF  YHW 

40  REM  nNO  FITCH  FOR  THE  FIVE 

50  REM  HOLE  PRESSURE  PROBE 

60  REM 

70  REM  DRV  ID  EVhNS    NOV  1586 

SO  RhM 

SO  REM  '...'PR  I  ABLE  HRMES 

i  6  0  R  E  M 

ilO  REM  E';:i>  IS  THE  UOLTRGE 

120  REM   RERD  FROM  THE  DRTR 

130  REM   RCQUISTION  SYSTEM, 

14  0  Re M 

150  RE M  fl  ■::  I  >  I  S  THE  U 0 N '•..' E R S  I  0 U 

ICO  REM  FROM  VOLTRGE  TO  INCHES 
155  REM   OF  WRTER. 

170  REM 

ISO  REM  P<:r:'  WHERE  1  =  1  TO  5 

lyO  REM   RRE  THE  PRESSURES 

200  REM   FROM  THE  5  HOLES 

210  REM   OF  THE  PRESSURE  PROBE 

220  REM 

230  REM  Pt"  =  STRTIC  PRESSURE 

240  REM  P7=T0TRL  PRESSURE 

250  REM  P:  S>=H'...'ERflGE  PRESSURE 

260  REM  Dl=  P2-P3 

270  REM  D2=P4-P5 

230  REM  D3  =  Pl-P':.i5> 

290  REM  Cl=  YRW  COEFFICIENT 

300  REM  C2  =  PITCH  COEFFICIENT 

310  REM  C3=CUEFFICIENT  OF  TOTRL 

320  REM      PRESSURE 

330  REM  C4=C0EFFICIENT  OF  STRTIC 

PRESSURE 

340  REM 

350  DIM  Cl<  1O0>  ..  C2<  108>  ..  C3<  100; 

360  DIM  C4  <  1 00  >  ,  y.  1  <  1 00  >  ..  Y '.:  100  > 

370  REM 

380  DISP  ^' ENTER  DRTE.  CMMDOYY  >  •' 

330  INPUT  Ml 

400  DISP  "ENTER  TIME.<HHMM>" 

410  INPUT  M2 

420  REM 

430  PRINT  Ml..  .■M2 

440  REM 

450  CRERTE  '•CRLDR1".2 

460  RSSIGN^H  1  TO  "CRLDRl" 

470  !  PRINT#  1  .:  M1..M2 

4  30  REM 
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4'^0 

500 

510 

520 

570 

540 

S'=i!.-1 

- 1-,  f1 


570 
530 
590 
f.  0  fi 


610 
620 

640 

650 
060 
670 

6  0  0 
690 

700 

710 
720 

7  3  0 
740 
750 
76^ 
770 
730 
7  9  0 

300 


310 

820 

330 

3  4  U 

350 

360 

870 

830 

390 

9  0  0 

910 

920 

REM  :HIS  3ECTI0N  COfipiirES 
REM  THE  C0H'...'ER5inN  FHrinp 
RtM  FOR  VOLTflGE  TO  INCHES 
REM  OF  i-JHTER.THE  USER  rili'^:T 
REM  INPUT  THE  MflHOriETER 
REH 

REM  THIS  SECTION  COMPUTES 
615  REM   THE  ZERO  OFFSET 'mp 
bit.  REM   RECTI  ON  FOR  THE 
TRHNS- 

REM   THE  ZERO  OFFSET  COR- 

622  REM   RECTI  OH  FOR  THE  TRfl 

HS-  623  REM 

REM 

DISP  "COMPUTE  ZERO  OFFSET  " 

DISP  "CORRECTION  .  •• 

DISP  "  ! ! !  DISCONNECT  fiLL 

TUBING  TO  THE  TRR 
NS  DUCERS  ill'' 

DISP  "  HIT  CONT.  WHEN  REfiDY" 
PhUSE 
CLERR 

FOR  1=200  TO  204 
Z  2  =  0 

FOR  .1=1  TO  50 
OUTPUT  709  .;  "Hi  "  .;  I  .;  "  '..'T  1  " 
ENTER  709  ;  X2 

NEXT^j"'^ 
Z3  =  Z2.''50 
Z4';  I-199>=Z3 
NEXT  I 
REM 
REN 

REM  :|:1c*:*::f.:*:*L00P#  1  *  :|c  .t  *  *  :*c  4:  *  :^ 
REM 

FOR  1=200  TO  204 
11=1-199 

DISP  "TRANSDUCER  CflLIBRflTION 

ENTER  THE  MflNOMETE 

R  PRESSURE  ':;iN 

H20>" 
DISP  "ENTER  THE  PRESSURE" 
D  I  SF  "  TRHNSDOUCER  NR  ."..II 
INPUT  H2 

Z0  =  0 

REM 

REM  t:^tttL00Pi/^2ttttttt.:if.t^: 

FOR  J=l  TO  50 

OUTPUT  709  .;  "HI  "  ..  I ;  " y T  1  " 

ENTER  709  ;  X 

Z0=Z0+X 

NEXT  .J 

REM  ttt-tf^rttttttttttttttt 
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"sfl 


r'  t.  y 

579 
9  8  H 
556 
1  Q  tj  9 
1  y  1  9 
i  9  2  9 
i  9  3  9 
1  9  4  9 
1959 
1  9  6  9 
1979 


Ei=Z9.-59 

H^I-15y>=H2^--<El-Z4<I-155V;. 
PRINT  "R"..  I..  "  =  "  ..  a<I-i55  J 
DISP  "IF  H(I>  IS  UHShT 

PRESS  S  " 
INPUT  H2 
IF  h2=S  THEN  GOTO  399 

he;-;t  I 

REM 

R  E  M  *  *  ;*;  1 1  :*:  i:  :*  t  *  *  *  1 1 1 1  %  %  %  %  X  %  X 
REM 

REN  ENTER  THE  'v'RLUES  OF 
REM   STATIC  HHO  TOTRL 
REM    PRESSURE . 
REN 

DISP  "INPUT  STATIC  PRESSURE 
< INCHES  OF  H2U>" 


1 9  S  9 
1999 
1  199 
1119 
1  129 
1  139 


1  1  b9 
1  179 
i  1S9 
1  1  59 
1299 
1219 
1229 
1239 
1249 
1259 


CIS 
DIS 
REN 
D I  S 
D I S 
OIS 
PAU 
CLE 
B6  = 
FOR 
REM 
OUT 
ENT 
B6  = 
HEX 
Ph  = 
REM 
DIS 


P  "USE  TRANSDUCER  NR  1 " 
P  "AND  CHANNEL  HR .  299" 

P  "HIT  CONTINUE  WHEN" 
P  "P  STATIC  PROPERLY" 
P  "HOOKED  UP. " 
SE 
RR 
9 
.J=l  TO  59 

PUT  795  ;  "  A I  "  .;  299 ;  "  '...'T 1  " 

ER  795  ;  K6 

a6  +  X6 

T  J 

Be-'-59 

P  "INPUT  THE  TOTRL  PRES- 
SURE. < INCHES  h; 


1260 
1279 
1239 
1259 
1399 
1 3 1 9 
1329 
1339 
1349 
1359 
1  369 
1379 


DIS 
DIS 
DIS 
DIS 
OIS 
PflU 
CLE 
REM 
B7  = 
FOR 
OUT 
ENT 


AND  CHRNNEL  NR . 299 

II 

HIT  CONTINUE  WHEN" 
P  TOTRL  PROPERLY^' 
HOOKED  UP . " 


P 

P 

P 

P 

P 

SE 

AR 

9 

J=l  TO  50 
PUT  795  .;  "RI" 
ER  795  .  X"7 


:90.:  "UTl  " 
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1339  B7=B7+X7 
139^J    MEXT  J 
1409  P7  =  B7.-'59 
1419  REM 
1429  REM 
1439  REM  ENTER 


RMBIEHT    CiJNDrTI9H 


14 
14 
1  4 
14 
14 
14 
15 

-■  c 

i.  J 

15 
1  5 


15 
15 
15 
16 
16 
16 
16 
16 


49  REN 

59  DISP  "ENTER  PRMBUH  .  HG)  " 

69  INPUT  P.l 

79  REM 

89  T=9 

99  FOR  J=l  TO  59 

99  OUTPUT  799  .  "  fi  I  "  .;  1 99  ;  "  UT I  " 

19  ENTER  799  .:  Tl 

29  T=T+T1 

39  NEXT  J 

49  T2  =  T.-'50 

59  GO SUB  3919 

69  REM 

79  REM  DEL  P  COMPUTED  FROM 

39  REM   PTOTaL-PSTRTIC 

99  REM 

99  r:  =  P7-P6 

19  REM 

29  REN 

39  REM  CONVERSION  TO  31  UNITS 

49  Hl=fll.t33o5  .  32 


io59 
1669 
1679 
1639 
1699 
1709 
1719 
1729 
1739 
1740 
1759 
1769 
1779 
1739 
1799 
1390 
1319 

1329 
1339 
1349 
1359 


C  =  C:t: 

T4  =  T 

R1=H 

U1  =  C 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

FOR 

K1=K 

REM 

REM 

REM 

DISP 

INPU 
REM 
DISP 
DISP 


243  .  7 
3+273 .15 
i/<2S7:t:T4> 
2*C.-R1>-V5 


ENTER  THE  BIG  LOOP 
FOR  RERDING  PRESSURES 
RND  COMPUTING  CuEFF . 

K"=l  TO  199 

ENTER  THE  VRLUE  OF  THE 
PITCH  RNGLE  RND  YRW . 

'•  ENTER  THE  VRLUE  OF 

THE  PITCH  RNGLE" 

T  X  1  V  K  > 

"ENTER  THE  '...'RLUE  OF" 
"THE  YRW  RNGLE  .  " 
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1860  INPUT  YvK::- 

1:570  DISP  "  PITCH  hHGLE  ENTERED 

I S  "  ..  X 1  (  K  ::• 

1530    DISP     "YHW    HHGLE    ENTERED     IS 

"  ..  V  <  K  > 

1-590    DISP    " 

1500  DISP  "  DO  YOU  WhHT  TO  CHHHG 
E  PITCH  OR  YRW  IF 

SO  PRE  SS  1" 

1310  INPUT  N4 

1520  IP  N4=l  THEN  1SO0 

1930  REM  THIS  LOOP  HCQUIRES  EhCH 

1940  REM  PRESSURE  50  TIMES  AND 

1950  REM  ThkES  THE  R'-^ERflGE . 

1  9  €•  0  R  E  M 

1970  DISP  "HIT  'CONTINUE'  FOR" 

19S0  DISP  "DflTR  fiCQUISITION . " 

1990  PRUSE 

2000  CLEhR 

2010  REM 

2020  REM  :t:t.;*:ic:1clcL00P#3-t:*1c:«c***:^;t: 

2030  FuR  1=200  TO  204 

2  040  Z 1 =  0 

2050  REM  *:4:i:.l::*;.tLOOP#4.*:**:t:*:1c:4:** 

2060  FOR  J=l  TO  50 

2070  OUTPUT  709  ;  "R  I  "  .:  I .;  "UT  1  " 

20S0  ENTER  70?  ..  X 

2090  Zl=21-t-X 

2100  NEXT  J 

2110  REM  4: ;<c.*:.t :*::♦: .*::4:1c:1c **lc.-*: :|r *:♦:** :t:*:t; :*: 

2120  X  =  Zi'-50 

2130  E<I-199>=X 

2140  P  <  I  - 1  99  >  =R  •::  I  - 1 99  ;•  .*:E  •:;  I  -  1  99  )  - 

Z4< I-199> 
2150  NEXT  I 

2160  RE  M  %  t  %  1 1 1  *  *  *  :*:  t  *  1 1 1 1 1 1 1 1  *  1 1 
2170  REM 

2130  p<ey=Q 


2190  FOR  1=2  TO  5 

2  "'00  P-:  6>=p<  I  j'+P-re) 

2210  NEXT  I 

2220  P<6>=P<6>.'"4 

2230  D1=P':;2;'-P'.:3) 

2240  D2  =  P':4:)-P<15> 

2250  D3  =  P(:  1  >-P<6> 

2260  C1(K>=D1/D3 

2270  C2CK;:'=D2^D3 

22S0  C3<K>  =  <PU  >-P7>.-D3 

2290  C4':K;'  =  ';P<6>-P6>''D3 
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2300 
2310 

2320 
2330 
2340 

2350 
2  3  €•  0 
2370 
23S0 
2390 
2  4  0  M 
2410 
2420 
2430 
2440 
2450 
2460 
2470 
2  4  80 
2490 

2500 
2510 
2520 
2539 
2540 

2550 
2560 
2570 


25S0 
2590 
2600 
2610 
2620 
2630 
2640 
2650 
2660 
2670 
26S0 
2690 
2700 
2710 
2720 
2730 
2740 
2750 
2760 
2770 


REM 

PRIHT#  1  :  ;:<1<:K>.  v<K>..  Clt;K> 

..  C2'- K>  ..  C3<K>  ..  C4';k  > 

REN 

RErl 

PRINT  ''PITCH  ANGLE  IS'-.^KK 

F  R  I  h  T  "  V  fi  W  R  N  G  L  E  IS"  ..  Y  <  K  > 
PRINT 

PRINT  "Pl  =  "  ..P<  IM- 
PRINT "P2="..  P<2> 
PRINT  ='P3="..  P<3> 
PRINT  "P4="..  P<4:j 
PRINT  "P5="..  P<5> 
PRINT 
PRINT 

DISP  "LRST  POINT'^?" 
0I3P  "IF  SO  PRESS  1" 
INPUT  Nl 

IF  Nl=l  THEN  GOTO  2500 
NEXT  K 

R  E  M  t  %  %  X  %  %  t.  %  %  t  X  X  %  X  X  X  X  X  X  X  X  X  X 
1705 
REM 

HSSIGN#  1  TO  X 
REM 
PRINT 

PRINT  "X X X X X X X X XXxX X X X X X X X X 
XXXXXXXX" 
PRINT 
PRINT 
PRINT  "FIVE  HOLE  PRESSURE 

PROBE  CflLIBRRT 
ION.  " 
PRINT 
PRINT 

PRINT  "DRTE  OF  RUN  IS"..  Ml 
PRINT  "TIME  OF  RUN  IS"..  N2 
PRINT 
PRINT 

PRINT  "DENSITY';K;G/M3>" 
PRINT  Rl 

PRINT  "'..'EL0CITY<M-'3>  " 
PRINT  Ul 

PRINT  "PRMBc:N/M2>" 
PRINT  Hi 
PRINT  "TRMB(C>" 
PRINT  T3 

PRINT  "P  STRTIC  -^IN.  H20>" 
PRINT  P6 

PRINT  "P  TOTRL=<;iN.  H20  .- " 
PRINT  P7 
PRINT 

PRINT    "PIT     {H\-\         CPl  c 

P2" 
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273H  REM  :t:**:*v^::^L0uP#5;t::*::*:*t;t:**:*;* 

2750  POR  1=1  TO  Kl 

2S^e    P-RIHT  USING  2310  ;  X  1  •:  I  >  ..  Y  < 

1  )  ■  C  i  '•!  I !-'  .'02  (  I  } 
23  1  e  i  MflGE  MOD ..  2:-: ..  MOD  ..  2X ..  MD  .  DDDD 

.2         X..MD.DDDD 
2820  NEXT  I 

2  330  REN  .^.t.^:^:rtt^ti:tt^-^ttttt*1lit 
23  40  PRINT 
2S"?0  PRINT 
2360  PRINT  "PIT   YftW    CP3     CP 

2370  REM  ;*.:***:t:tL00P#b*:*:******4;:ic 

2330  FOR  1  =  1  TU  1(1 

2S90  PR  INT  US  I NG  2909  ;  >=!  <  I  >  .  Y  < 

I  :■  ..  US-;  I  >  ..  C4v  I  > 
2300  I  MfiGE  MOD  ..  2X .  NDD  ..  2X ..  MD  .  DDDD 

.2         X.ND.DDDD 
2910  NEXT  I 
2920  REM  :t:*:*;|::^-:*::*:t:t:t.t**. let-* :*:*:* ***:!: 

:*: 
2930  PRINT 

2940  DISP  "  THRNKS  FOR  USING  ME" 
2950  DISP  "  HOPE  YOU  FOUND  WHftT" 
2960  DISP  "YOU  WERE  LOOKING  FOR" 
2970  DISP  "  REGRRDS" 

2930  DISP  "  HP-85e" 

2990  END 
3000  REM 

3010  REM  yOLTRGE  TO  TEMPERRTURE 
3020  REM  C0H=.;ERSI0N 
3030  REM 

3040  REM  yfiJVJy     TO  T-lO 
3050  REM 
3060  E1=T2*1000 
3070  T3  =  26.573*E1-1  .  936879:?iEl  .*E1 

+ .99785*E1*E1*E1- .2bl277*El 

:t:El*EitEl 
3030  RETURN 


SI 


20 

REM 

30 

REii 

40 

REM 

50 

REM 

60 

RE!^! 

70 

REM 

88 

REM 

30 

R£n 

1  'i  0 

REn 

i  LG 

REh 

IJQ 

F;£;i 

L3S 

REri 

i4G 

REM 

150 

REM 

T  l"  •"■! 

REM 

170 

RE!^, 

t  i' r- 

R  R  f'i 

1S0 

REM 

ZOQ 

REM 

Z3i5 

R£ri 

244? 

REh 

-,  c?  b 

REM 

260 

REh 

270 

REH 

"'  P  f"" 

n-B'. 

292 

REM 

300 

REfi 

pCM 

320 

REri 

REM 

340 

REri 

353 

REN 

360 

REM 

370 

REM 

380 

REM 

330 

REM 

400 

REM 

410 

REM 

420 

REM 

430 

REM 

440 

REM 

450 

REM 

450 

REM 

470 

REM 

480 

REM 

4  5"0 

REM 

']00 

C'CM 

S10 

kF.M 

S20 

REM 

rn'J'--nr,n    r-Rsfv.yU 

THIS  RRCGRAM  nCQUIRES  THE  PRESSURES  FROM 
A  FIvE  MCi.E  RRE;SURE  PRUEE  AfJO  COMMUTE:: 
OBTAINED  FROM  TiiE  PRObRfiM  "RROCAL'' 

AUTHOR  DAVE  E^AN;^ ,  JArjlJAR  / ,  FESRUAR  Y  ,  13(^7 

I/ARIASlE  rJAMES 

E(  I  }  IS  THE  yOLTAGE  READ  f-RON  THE  DATA  ACQUISTIOrJ  £\STEM 
A(I)  13  THE  CONUERSION  RROil  UOLTnEE  TO  ITjCHES  GF  WATER 
R(  I  )  13  THE  PRESSURE  READ  FROM  Ti-IF  fi\!f    HOLE  PRESSURE  PROBE 
1=1  TO  5  AND  IN  EACH  OF  THE  ABOVE  IS  ASSOCIATED  wITH  THE 
PRESSURES  FROM  THE  FIl'F  HOLE  PROBE 

PS-STATIC  PRESSURE  OF  THE  FREESTREAn 
R7=T0TAL  PRESSURE  OF  THE  FPEE5TREAM 
R  (  S  .'  -  A  V  E  R  A  G  E  P  R  E  S  S  U  R  E 

D2=F4-Pb 
D3=P1 -R(b ) 

Ci^GAUlBRATED  YAW  COEEFICIENT 
C2=CALieRATED  PITCH  COEFFICIENT 
C3=CALIBRATED  TOTAL  CQEFFICIENT 
C4==CALiERATED  STATIC  COEFFICIEriT 

MI-MEASURED  YAw  COEFFICIENT 
M2=MEASURED  PITCH  COEFFICIEMT 
M5-MEASURED  TOTAL  COEFFICIErjT 
M4=MEASURED  STATIC  COEFFICIEflT 

^-'=^THF  LuCAl  TOTAL  v^ELOCITy  UECTOR 
l'l=THE  X  COMPONENT  OF  VFLOCIT^■ 
U2-THE  Y  COMPONENT  OF  UELOCITY 
UZ=THE  Z  COMPONENT  OF  UELOCIT; 

Xl^THE  PITCH  ANGLE  FROM  THE  CAlISRATIuN'  OATA 
y-THE  YAW  ANSLE  FROM  THE  CALIBRATION  GATn 

X4=THE  PITCH  ANGLE  OF  THE  FLOlv 

Y1  =  THE  VERTICAL  LOCATrON   OF  THE  •^RCBE    TIF 
Z0=THE  SPANWISE  LOCATION  OF  THE  PROOF  TIF 
Y3-THE  YAW  ANGLE  OF  THE  FlQU' 

ril -AMBIENT  PRESSURE 
T4-ABS0LUTE  AMBIENT  TEmP( KELVIN  ) 
RI=DENSITY  (KG/M-3) 
UB-UElDCITY  of  THE  FREESTc■r:?i^'| 
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530 

540 

570 
eg '3 

593 

£00 

bie 

620 
S30 
S40 
650 

6'^  3 
680 
S3S 
70? 


740 
750 
750 


813 

320 

840 

S5C 
o  c  '■> 

Q7i7i 

680 
890 

gi0 

920 

93Q 

340 

950 

36  0 

3d  1 

370 

980 

323 

1300 

1018 


DIM  Y(  1000  ),Z(  1030  \t^i\   1000  ^n2(  1000; 
DIM  Pl( 1000), PB< 1000  ; 

REr-1 

I^;^•■JT  "EMTER  DnTc: ,  (  ^'^■DUV•/ )  %r^8 

INPUT  -hNlhR  Tln[\  ;HHMM  ;%n;i 

PR IN1ER  IS  701 

PRINT    "DnTE    OF    RUN    IS" , U8 

PRIfIT    "11  TIE    OF    RUN    IS"'  ,rJS 

KcPl 


[NKiir    "ENTLR 


INPUT  -Er-nER  PT5  '.'ERT  ICAL"  ,  N3 

REH  N3  nuST  BE  hii   EUEN  I^JTESER 

INPUT  "ENTER  SPnMUilSE  RESOLUTiGN(  IN  ; "  ,  24 

I'^PUT  "ENTER  VERTICPiL  RESOLUTION^  IN  )  •' ,  V4 

INPUT  -INITInL  2(  INr'  ,Z3 

INPUT  '•  INITIAL  v;  IN  :"  ,  y3 

PRINT  "PIS  SPANUISE"  ,r!3 

PRINT  "PTS  VERTICAL" ,N3 

PRINT  "SPANUilSE  RESOLUTION  ;IN;",24 

PRINT  -UERTICriL  RESOLUTION  (Ifj)",y4 

PRINT  "INITIAL  Z( IN)" ,23 

PRINT  "INITIHL  Y(  IN}-',Y3 

REN 

Zi 1  )  =  Z3 

Y' 1  !  =  Y3 

N9=N3/2 

PCR    17=1    TO    N9 

I G  =  1 7  -  i 

Jl=2Mb*n3r2 

J2  =  2'^IS*M3  +  ri3 

POP   K=Ji    TO    J2 

2Cr  >=Z>  K-1  >  +  Z4 

Y' K  ;  =  Y(K-1  ) 

NEXT    K 

J3=2*IG»H3+N3+1 

2( J3  i  =  Z( J3-1  ) 

Y' J3;  =  Y<  J3-i  )  +  Y4 

J4=J3rl 

J5=2»I6*M3+2*M3 

FOR  K=J4  TO  J5 

2(K  •  =  Z'K-1  )-Z4 

y  ( !<  )  =Y  (K-1  ) 

NEXT  K 

IF  r7=N3  THEN  iQ10 

JS=J5+1 

Z( Jb  )  =  Z( J5-1  ) 

r( JB  )  =  Y( JS-l  H  v4 

NEXT  17 

REN 
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JRF^CnON 


10^0      ASSIGN   §Path2    TO    "PRSZ" 
1031       INPUT    "DO    YOU    WANT    TO    SKIP    FREi^lnirlnP. 
1332       IF    N$--="Y"    THEN    ZZS'Ti 
1040      REn 

1050      REn    THIS    oECTJOM    CQMp!jTE5    7ri£    ZE-0    Of 
ieb0      REn    FOR    THE    TR-iNSDUCERS . 
1070      REM 

1C1S0      FRIMIER    lb    1 

1050      pwINf    "i'-OnPUTE    THE    ZERO    OFFSET    CORRECTION" 
1100      PRINT    ■■  1  !  i    OISCO'JNECT    -^LL    fUBiNG    10    THE    TPriNSDLiCEP; 
1110      PRINiT    '■    HIT    CONTIMIJE    UHEr-i    PEPiDY" 
1120      PAUSE 

1130      fOf(    1  =  200    TO    204 
2-0 
OR    .1-^1.    TO    E0 

.UTPUT    709;  ''nr'  i  I  ;  H^Tl'' 

EN  i  iiR    TicjS  ;  X2 

G2^G2+a2 


i  (  ,  ■■ 


1140      f 
1150       r 

1150     : 
1170 
1180 
1  i  30 
1200 

1 2  1 0 


'.t?! 


1240 

1250 

»  *">  r*  /■> 
i  i_  C  -O 


1280 
12S0 

1300 
13i3 
1320 
I  u  J  •a 

1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1420 
1430 
1431 
1432 
1440 
1450 
1460 
14  70 
'4'd0 


PRESSORE    If!    < 


04(  I-lOg  ;  =  G3 
NEkT    I 

PEri    End    OF    LOOP   FOR    ZERO    CORRECT lOM    OF    TRriNSDuOERS 
REM 

REn    THIS    LOOP   O^LIBROTES    THE    TF<nN5DO0ER    rtSnINST    THE    'lANOnETER 

RErl 

FOR    1=200    TO    204 

1 1  =  1-159 

f^PINT    ■■    TReipiSDOOER    CfiLIBRriTION,    ErjTER    THE   MriNOnETEP 

U.ATER" 
PRINT 
PRINT 

PRINT  •■ErjTER  THE  PRESSURE  FOR  rRANSDUOER  NR .  ",1  i 
INPUT  H2 
G0  =  0 

FOR  j=i  TO  50 
OUTPUT  709;  '^nl"  ,  I;  "OTl" 
ENTER  709; X 
G0=S0+X 
NEXT  J 
Ei=G0/50 

ni  1-199  )=H2/(E1-S4( 1-199  ;  ) 
PRINTER  IS  701 
PRINT  "AM,  ■'  =  '•, r!(  1-199) 
PRINTER  IS  1 

PRINT  "A" , 1, "=" ,^{ I-1S9 ) 

INPUT    "IS    A( 1-199)    SHTISFnCTORY    (V/NJ-.NS 
IF    NS=-N"    THEN    1230 
PRINT 
PRINT 
NL''T    I 


'4CHES    Of 
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1500 
1510 
1520 
1530 
1543 
1550 
15B0 
i57G 
1530 
1530 

1S00 

1610 
1S20 
1630 
lb4e' 
IB  53 
1660 

i  b  ''0 
.  -->  #-k  r-» 

1590 
1700 
1718 
1720 
1730 
1740 
i  753 
1733 
1770 
1730 
1790 
1808 
iSiC 
1820 
1830 
1340 
1S50 
1360 
1S70 

i    .~,  .-V  /-t 

13S0 
1300 
1910 
1920 
1930 
1940 
1350 
1958 
1370 
1330 
1990 
2000 


REn    END   OF    LOOP    FOR    CALIBRATItNG    TR^•M5iJUC■EPS 

REM 

REn   ENTER    THE    v'ALiJES    OF    ST/iTlL    nNO    TGTnL    PRESSURE    OF    THE    FREESTREf^M 

PRINT    "INPUT    STnTIC    PRESrURE    OF    FREEi^TREAM    ^IN.H20)" 

PRINT    "USE    TRANSDUCER   ^iK .  1    rViiD    CHr^NnEL    rjR.200- 

PRINT 

PRINT  "HIT  CONTINUE  WHEN  P  hTriTIC  IE  PROPERLV  SET  UP" 

PnUSE 


■  L^r, 


I  u  b 


C?,'! 


CLiTPUT  739;  "AI"  ;209;  ••UTl'- 

ENTER  709; KB 

BS=6S+X6 

NEXT  J 

REM  END  LOOP  FOR  ACQUIRING  3TnTIC  PRE^2. OP  FREEiTRErirl 

P9='36/53 

PRINT 

PRINT 

REl^' 

PRINT  "INPUT  TCTr^L  PRESS'JRE  OF  "hE  -REESTREriM  ,  (  IN  .  h20  ) " 

PRINT  "0!£E  TRANSDUCER  NR .  1  AND  CHANNEL  NR.200" 

PRINT  -^MIT  CONTINUE  WHEN  P  TOTAL  PROPERLY  SETUP" 

PAUSE 

B7=0 

FOR  J=l  TO   50 

OUTPUT  70S;  "AI"  ;  2Q8  ;  "UTi  •• 

ENTER  70S;X7 

B7=B?+X7 

NEXT  J 

REfl  END  OF  LOOP  FOR  ACOUIRIr^G  TOTAL  PRESSURE  OF  FREESTREAN 

P7=B7/50 

REM 

REN  ENTER  ANEIENT  CONDITIONS 

INPUT  "ENTER  PANB  iIN.OF  Ho!",Al 

RErl 

REM  LOOP  FOR  ACQUIRING  Tb.nP    OF  FRFESTREAN  >JIA  rHERNOCOUFLE 

1  =  0 

FOR  J=l  TO  50 

OUTPUT  70S;  "AT'i  109;  "UTl" 

ENTER  703; Tl 

T=T+T1 

NEXT  J 

REM  END  OF  LOOP 

T2=T/50 

REM  CONVERT  uOLTAGE  TO  TE^iPEFATUPE 

Ei=T2*1000 

T3  =  2d.573»E1-1  .  335o7S'^El  •■  2+ .  33785*E  1  "5- .  251 277  ♦El"  4 

REM 

REN  CONUERSION  TO  SI  UNITS 

T4=T3+273. 15 

nl=Al *  3385. 82 
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2010 
2320 
2030 
2040 
2050 
2060 
2070 

2330 
21iSC 
2  113 
212ii 
2130 

2  L4.0 

2150 
2ib0 
2170 
2  180 

-~  1  rv  r\ 

2203 

'-  -  i  -s 


2240 
2250 
2260 
"^  '•'  ~  1/1 

2280 

2300 
2331 
2  jk'i- 


R1=^1/C28?*T4 ! 

C  =  (  P  7  -  P  9  )  *  2  i  S  .  7 

REH  FREESTRErifl  VElOCITY 

U9  =  ',  2*C/R1  >'^  .5 

REfi 

REM  BEFORE  £P\C'r^    Dr\Tf\    RUN  ENSURE  THnT 

PRINT  "EnLANCE  P2  AND  P3.  " 


THE  PROBE  IS  BnLANCLO 


2330 
2331 
2340 
2390 


2  4^y 
2430 
2440 
24S0 
2468 
2470 
2430 
2490 
2530 


PRINT  "HIT  CONTINUE  WHEN  ALL  TuBUio  HnS  GEEN  CONNECTED" 

Pr^USE 

REM  LOOP  FOR  BALriNCINS  PROBE 
FOR  1-231  TO  202 

FOR  J=:  TO  50 

OUTPUT  703;  ''AI  "  ;  T  ;  "UTl  " 

ENTFP  709; K 

G  i  =G1  +'^' 

NEXT  J 

X=Gi/50 

El  1-139  )~)i 

Rv  I-iSy  )  =  rta-  139  )*E<  I-I33  ;-G4(  1-199  ) 

NEXT  I 

PRINT 

P  R I  iN  T 

PRINT  '■P2=vp{2) 

^RINT  ''PO^"  ,  r  (  3  ) 

INPUT  ^•GOES  P2=P3  (Y/N}',N!£ 

IF    M$="N"  THEN  2100 

PEN 

REN  ENTER  THE  LOOP  FOR  f=)CQUIRiNG  EnCH  PPESSuRE  OOiiPuriNG  CGE:i^FlCEinS 

REN  AND  COMPUTING  TUQ  PRESSURE  COEFFICIENT:  nNO  TWO  PPE32URE- 

PRINTER  IS  701 

PRINT  -y     Z      Ml      h2  PI      Pril^G" 

PRINTER  IS  1 

REN 

K9  =  a 

i;2=N3>N3 

FOR  K=i  TO  K2 

K9==K 

REN 

WnIT  10 

REM  LOOP  FOR  ACQUIRING  EACH  PRESSURE 

FOR  1-230  TO  204 

61=0 

REM 

FOR  J  =  l  TO  10     !  ENTER  THE  DAS  AND  SAMPLE  EACH  PRESSURE  tO  "^IMES 

OUTPUT  709;  "AI"  ;  I  ;  "V'Ti" 

ENTER  7S9;X 

Gl-oi+X 

NEXT  J 

REM 

X-Gi/10  lAUERAGP  THE  VALUES 
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'';:'t7i 


2553 

:sB0 


E(  I-19S '  =  A 

F(  1-199  ;=f^<  1-133 '*E;  1-139  ;-G4',  I-lG-3  ■     iCOPPtCTION    FOP    ZERO    AND    TR/=,rjS,    CrI 

IF  i-:30  the:] 

PnK  )-P(  1-199  > 
END    IF 
NEXT    I 
PEr^   END    OF    LOOP   FOR    -^CCUiRIrjG    ^RES5U'-E^ 


REn  CUMPUTE  THE  nvERnOE 


.■^c 


t-^.- J,^-i,^•■- 


2530 
2530 
2530 
2  S  1 0 
<_  b  ^  ii/ 
2S.30 
2  =  40 
2E50 
2551 

'"» f"^  r^  -"i 


2710 

'J  '^  "^  iT^ 

2^40 
27E0 
2760 
2770 
2780 
273Q 
2980 
2  3  i  0 
2320 
2830 
2340 
2353 
23G0 
287Q 
2339 

-1  -■  n  n> 
i-O  J<cJ 

2308 
2310 
2323 
2330 
2340 
2350 
2  3S0 


1J6=P(  I  .'-^05 

NEaT  i 

Po'K  )  =  0G /'i 

-EM  END  LOOP 

r-;Ei1 

REi"!  CGMFl'TE  THE  COEFFICIENTS  Or  VnuJ  n'-iD    PITCH 

RErl 

01-F( 2  )-F( Z  ■ 

D2  =  Pv  4  )-P(5  ) 

C3=F( 1  ;-PG' K  ) 

Mi(K)=Dl^'05       i  LOCAL  V'r^vJ  CCEFFICETiT 

r!2'K;=D2/03      'LOCAL  PITCH  COEFFICEraT 

FPINTEP  IS  701 

PRINT  USING  2750;  Y\  K  ; ,  2;  K  /  ,rll- K  ;  ,i^2(  K  ^  ,P1:K  )  ,PS\  K  ; 

lr-!AGE  MQO .  CD  .  2a  ,  nCD  .  DD ,  LX  ,  riDu  .  DDDD  ,  2X  ,  MOD  .  DOCD  .  2  >: ,  nDD  .  ODDD  ,  2X  ,  flOD  .  DGDD 

P£ri 

NEXT  K 

REM 

FOR  1=1  ~0  K9 

OUTFUr  5P5th2;  Y(  I  ),2'  I  ;.i^i:  I  ).'':2(  I  )  ,F  I  ^  I  ;  ,  PS'  I  / 

NEXT  I 

H33IGN    ©Path!    '^0    * 

.ASSISr-J   @Pa th2    TO    * 

Pj^'jfjj     '»•♦<<»♦**■«■>*»♦»  »»«-»*)i**»»-i- •■«••■♦+*«-•■»-»♦»»»♦••♦»♦**  +  «♦♦♦" 

^•RINT 

PRINT  "VELOCITY    COflPUTAT ION   Ar-:C    v'ELOCITY    COrPOf^ENIG  ■' 

PRINT 

FRIiiT 

PRINT  "DATE    OF   P.UN    I5",ri3 

PRINT  "TIME    OF   RUN    IS'',N9 

PRINT 

PRINT  " DErJS  I T Y  i  K Q / M -  3  ) "  ,  H 1 

PRINT  -FPEESTREAn   UELOCITv    (rl-3;",U3 

PRINT  -PAMBIENT(N/M-2  :- .Al 

PRINT  '■TAnBIENT<C)%T3 

END 
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20 

4  0 

50 
60 


33 

lis 

120 
130 

U0 

150 

iB0 

173 
180 
130 
200 
Zli? 
220 
230 
240 
250 

n  c  .'^ 

L.  U  iy 

270 
280 
?9P 
300 
310 
320 
330 
343 
35G 
3S0 
370 
38Q 
390 
400 
410 
420 
430 
440 
450 
450 
470 
480 
4  30 
S00 


REn  PRCGRnM  vEL 

Ren  THIS  PROGRnrl  COr-IFUTCS  THE  COMPONENTS  OF  VELOCITY  FROn 

RErl  THE  RnW  DATA  COLLECTED  v'lA  THE  PROGRAM  PRSACQ 

REM 

REM  AUTHOR  uAv^'E  EVANS,  J  ANUARv  ,  FEBRUi^R  V  ,  1987 

REM 

REM    I'ARIABLE    NAMES 

REM    P5=AyERA3E    PRESSURE 

REM 

REM  rii-MEASURED  VAU  COEFFICIENT 

REM  M2=MEn3URED  PITCH  COEFFICIENT 

REM  M3=MEn3ijRE0  TOTAL  COEFFICIENT 

REM  M4  =  MEASL!RED  STATIC  COEFFICIENT 

REM 

REM  U=THE  lOrAL  VELOCITY  AT  THE  pi^'OBE  TIP 

REM  Vi-TH£  K  COMPONENT  OF  iJELOCITY 

REM  V2-THE  V  COMPONENT  {^.'r    UELOCITY 

REM  vZ  =  THE  Z  COMPONENT  OF  I'EtOCITY 

REM 

REM  Xi-THE  PITCH  ANGLE  FROM  THE  CALIBRATION  DATA 

REM  Y-THE  YAUi  ANGLE  FROM  THE  CALIBRATION  ri.-.Tn 

REM 

REM  >'.4-THt  PITCH  ANGLE  OF  THE  FLOW 

REM 

REM  Y  =  THE  VERTICAL  LOCA'aON   OF  THE  PROPE  TIP 

REM  Z-THE  SPANWISE  LOCATION  OF  THE  PROBE  TIP 

REM  Y3=THE  yaw  ANGLE  OF  THE  FLOW 

REM 

OIM  Ci*  77  )  ,C2(  77  )  ,C3w7  i  ,C4'  77  )  ,Mi(  830  ).M2;  000  '< 

DIM  Xl(77),  Y{77  },Y  1(800  ) ,  P  1  (  SO0  ) ,  PS(  800  S  PS '300 

DIM  Z0f  800  ),  11(800  ),\;nS00  )  ,ij2':S00  ),'.'3(8S0  ; 

DIM  2^800  ■,P(S00  ),U4(e00  ),vb(300  !,VB^300  Sv'Z^SOO) 

REM 

INPUT    "WHfiT    13    THE    DENSITY    FOR    THIS    RUN'!  KG/rr'S  ) "  .  Rl 

REM 

REM  READ  CALIBRATION  DATA  INTO  COMPUTER  M£MOF\ 

REM 

ASSIGN  QPathl  TO  •'CAL" 

FOR  1=1  TO  77 

ENTER    apathliXK  I  ),Y>  r  ^Cl'  I  !.C2^  I}\03v  I  ?,C4Y  I; 

NEXT    I 

REM  READ  RnU  DATA  INTO  COMPUTER  MEMORY 

ASSIGN  'SPathZ  TO  ■■PR52" 

FOR  1=1  10  800 

ENTER  @Path2;  Yl(  I  ),Z'  n,Ml(  I  ),M2<  I  ),PU  I  ),PS(  I  ' 

NEXT  I 

REM 

CREATE  BOAT  "UEL2''  ,225 

ASSIGN  ^Pdtn3  TO  "UELC' 


S8 


c  ~  pi 

530 
540 

550 
563 
570 
530 
590 
bOQ 
die? 
5:3 
630 
640 

6S0 

670 

s  y  Q 

S90 

7,10 

r  i-W 

730 
74B 
750 
7S3 
770 
780 
799 

900 

810 
820 
S3G 
b40 
850 
SSiZ' 
873 
930 
893 
S08 
910 
920 
330 
340 
950 
960 
370 

yes 
3  3  a 

1Q03 

[.<J  i.  V 


FOR  K-I  TO  800 

VlvK  )-Vl{K  )*2.54 

2  '  K  /'  —  ""/'  K  '■  *  ^^  . ;;  4 

REn 

REp!  CunFuTE  THE  FIRST  APPROk  IHAT  IOn  OF  fAW  riPiGI.E 

B3=-. 15385924 

B4=-'^.35SG57I3       !  COEFFICENTS  nND  i^OLYMGrilpr-^^^L  i)RE 

C5-. 13550589        !  FROH  FORTRriM  PROGRAM  "POlV"  ANO 

CG-.30401553S        !  ARE  COMPUTED  FROM  THE  CP  '-nW  './S, 

C7=.0330000G9  !     YAU    DATri    FROM    THE    CriLIBRATION 

CB=- .330GS367 

y3=e3-S4*i-u vK  ;+C5»ni(K  )-2-^cs*ni; k  ;"3+c7*rn( k  .)-4-i-ca*Mi(K  )"S 

REM 

REil    GO    TO    TiiE    SUB-ROUTINE    FOR    PITCH    r;NGL£    IrjTERPOLnTION 

305UB    i40i? 

REM 

REM    GO    TO    THE    2UB-R0UTINE    FOR    VfiW    Ar-iSLE    INTEHFOLATION 

GOSUB    iG33 

REM  GO  TO  THE  SUB-RQUTlNb  FOR  OF  TOTriL  hNC  CP  STATIC  INTERPOLATION 
GOSUB  1350 

r^  c  ;*; 

REri  THIS  SECTiOM  COMPUTES  THE  MAGNITUDE  OF  THE  LOOAl  VELOCITY 

REM  VECTOR  AT  THE  PROBE  TIP 

REM 


Plih ^=P1(K  '*248.3     ICOMOERT  PI  TO  N/M '2 

PB(K )=PS(K ;*243.3     'CONUERT  PvS)  TO  N/M": 

REM  CALCULATE  STATIC  AND  TOTAL  PRESSURE 

P8=PS(K  )-M-i*;P;  <K  )-Pe(K  ;  > 

FR=-P8     IDONE  BECAUSE  INTERPOLATION  RESULTS  IN  NEGATI'JE  l;ALUE  FOR 

?vK  )=Pl{K  )-f'i3'vpUK  ;-PSCK  )  ) 

U4(K  ;-v2-^Ri>' PvK  ^-P3  )  )  •  .5 

REM 

REM  COMPUTE  THE  CQMPONEMTS  OF  UELOCITy 

LEb 

v5f  K  )='J4(K  )*C0S^X4  )*COS^  V3  ) 

^.'5'  K  )=U4(  K  )*SIi-!(  X4  ) 

;j7(K  >  =  U4<K  >*C0S':X4)*SIN(V3  ) 

NEXT  K 

REM  END  OF  DATA  REDUCTION  LOOP 

REM  RE-ORDER  DATA 

11  =  I 

12  =  40 

FOR  1=1  TO  72 i  STEP  80 

J  i  =  I  i 

J2  =  I2 

FOR  K-l  TO  43 
zo;  Jl  ;=2-  J2  ) 
P3^  Jl  ;=P( J2  : 
'J-'.  J  i  )--'.'^\  J 2  ■ 


PSTAl 
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1020     ui;  Ji  -^usv  j::  ) 

103i3  v^2(  j:  ;-UG(  J2  ) 

1040  U3( Jl  )  =  U7'  J2  ) 

1350  Ji=Jl+i 

1060  J2=J2-1 

1070  NEXT   K 

ii380  11  =  11+80 

1030  12=12+90 

1100  rjEXT    I 

1110  13-41 

1120  FOR  1-41  TO  8C3  STEP  Gf; 

1130  K=I3 

1140  FOR  J=l  TO  40 

1150  Z0(K  )-Z(K  ) 

11G8  pg(K/=F{K) 

ii7D  v(K  *  =  v'4(K  ) 

1  ISO  Ul  (  K  )-^Jb'_  K  } 

1130  v^2(K-=U6(IO 

1200  i-'SCK  i-U7(K  ) 

1210  K-K+i 

1223  NEXT  J 

1230  13=13+30 

1240  MEXT  I 

1250  FOR  K=l  TO  KS 

12c0  OUTPUT  §Pat 

1273  NEXT  K 

1280  r^SSIGN  gPathi  TO  * 

1290  ASSIGN  ©PathZ  TO  * 

IZ-iQ  ASSIGN  @Path3  TO  ^ 

1310  PRINTER  IS  I 

1320  PRINT  ••Y(Cr1)      7 

VZiM/S  )" 

1330  FOR    1=1    TO   K9 

1340  PRINT   USING    1350  ;  v  i  a  > .  ZOU  ; ,  PS^  I  ) /;•  1  ;  ,Vi ;  r  ^ .  m^.-  r  )    ^r^(i) 

1350      IflnGE   !1DD.DD,2X,nDD.DD.2A.M0nn   OOP !  2^' Vnn*  nnn  " -V  i^n;-,"  nnn  ' 
.ODD  ■  ■-       .  -    .      ^,-.>,  .......L... 

1363  ME>vT    I 

1370  PEM 

1330  STOP 

1390  REM 

1400      REM   PITCH   ANGLE    iNTERPOi riT lUN 

1410      REM 

1420      FOR    1=1    TO    77 

1430      IFY3<y(I)    THEN    1450         i  SEriRCH    THRU   CALlPP^iinri   n^^TA    r-,    ctm 
1440      NEXT    I  •■  -      u  .,,,.,  i,. 

1450  K1=I 

1460  K2=I-7 

1470  K3=I-i 

1480  !:4-IrB 

1430  FOR    j=K2    TO    k3 

l5iZi0  IF    h2(Hi<C2'Ji    THErJ    IFZZ 


t- 1  u  i  ( r^ 


I ; .-  M  ; 


UX(!-1/S 


'JVCM/S 


(,nDD.DDD.2 


u 


!VHW  A^iGLhS  WHICH  BRACKET  THE  r^PPROX  YAW 
i ANGLE 


'rlr-;0  THE  PITCH  HNGLE  WHICH  CORPESPONOS  TO 
t  THE  MEhSL'PEO  PITCH  COEFF. 
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1EI0 

15:0 

1533 

1540 

1550 
1560 
1573 
1580 
ISSi^ 

isii^ 
1 


bZ^     REr 


NEXT    J 

Zl  =  cc2(J  ;-r-i2(K  ))/*  C2(  J  )-C2(  J  -i  ; ) 

X2  =  X1':  J  i-21*(  Ai(  J  >-;<l^  J-i  )  )    'BRACKETING  PITCH  ANGLE 

FOR  j=r.  i  TO    K4 

IF  MZvK  •  :C2(  J  )  THEN  15  7ii 

f-iEXT  J 

Z2  =  ':  C2i  J  '-M^^  K  ■  )/;  C2C  J  :-C2'  J-I  ■  ' 

a3=ahJ  ;-Z2*(Xi(J  '-xicJ-i  }  ; 

2  3  =  ^  '/  (  K  4  ^  -  V  3  )  /  (  V  ■;  K  4  ;  -  Y  (  K  3  5  ) 

X4=X3-Z3*\K3-X2 ) 

RETURN 


!  BRACKETING  PITCH  i^NGLE 
i FITCH  ANGLE  Or  THE  FLOM 


15  38 
1540 

1S55! 
iS53 

Id  .  \L' 

1660 

16  SO 
1709 
1710 
1723 
1730 
174-0 
1750 
1753 
1773 
1 780 
1793 
1802 
1  a !  '"^ 

i  O  i  t/ 

1 3  2  0 
iS32 
1640 
1350 
1862 
1373 

isso 

1390 
1902 
1913 
1920 
1930 
1940 
1S59 
1960 
1970 
1960 


REf^  YAW  ANGLE  INTERPOLATION 

REM 

PQR  1=^  TO  77 

IF  V3^Y( I  -  THEM  IS80 

rJEXT  I 

f-j  "^  =  y  — ''.' 

H3=I-1 

PQR  j=H2  10  ri3 

IF   K4  :.xi:  J )  Then  I'^zid 


I 


iFINO'VAW"  nr^(=LE3  'AHICH 
iCORRESFONu  rO  THE  tlEiSSilf-Eu 
iCP  YAW  ON  PITCH  ANGLES  WHICH 
IE-RACKET  X4 


F3-J+G 

F4=J+7 

24-(  Ci  '  Fl  ;  -ni'  K  )  ;/ ^  CKFi  )-Ci^F4  )  ) 

y4  =  Y(F1  )-Z4*(  Y(  Fl  )-Y{F4  )  ) 

Z5  =  ^C1(F2)--H1<K  ;  //<Ci(F2  )-Cl(^3  )  ; 

V5  =  v(F2  )-Z5*(t(F2)-YsF3)) 

rS  =  (Xl(Fl  )-X4)/(XI(F1  )-X1(F2)-' 

Y3=V4-Zb*< Y4-y5 ) 

RETURN 

REM 

REn  INTERPOLATION  FOR  LOCAL  COEFFFICIENTS  Or    STATIC  AND  TO' 

Z3=0 

72  =  0 

FOR  I- 1  TO  77 

IF  y3'Y' I  )  THEN  1310 

NEXT  I 

F1  =  I 

F2=I-7 

F3=I-1 

F4=I+G 

FOR  J  =  !^2  TO  F3 

IF  X4;Xi( J  )  THEN  1980 

NEXT  J 

F5  =  J 

F6=j-i 

F7=Jr? 


■>tib  jC'i^c 
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2-210 
2QZ0 
Z030 
2  3  "f  kj 
2050 

2080 
203'3 

2100 

'^ '  1 1"^ 

L.  1.  i  k; 


Z3  =  (  Xl^FG  ;  -X4  /  •'(  Xi\  F6  ;-Xl(  FS  )  ) 
N4-C4nFd  ;-Z3*\C4; Fb  )-C4cFb  )  ? 
N3=C3; F6  i-Z3' ( C3( F 
L4  =  C-CF3  )- 


■  PF  ■ 


'_  -<  •  r  s 


L. -■  =  '. 


C3' F7) 


:3(ry 


12  =  !  V(  F6  ■->  3  )/■,  Y(  Fb  )-V;  F9  )  } 

ri4-rJ4-ZZ*'  N4-L4  }       !  LOCAL  STnl"IC  CCEFFI 

M3-N3-Z2*v N3-L3  )      i  LOCAL  TOThL  COEFFIC 

hfturn 

STOP 
£riD 


INf 
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1  REM  PROGRAM  PLOT 

2  REH  THIS  PROGRAM  r-^PiKEB  CONTOUR  PLOTS  FOR  STREnMVISE 

3  REM  '^ELOGITY  AArJD  CAN  BE  JlODIFIEO  10  MAKE  CONTOURS  PLOTS 

4  REM  FOR  TOTnL  PRESSURE 

10  Dlri  vii  8S50  ),Z0<  80Q;,P9(e00  ),v'303  ■  ,L1(90O),U2(  SSS  !,v'3(Sa-0 

20  DIM  X! SQ3  ) 

30  ASSIGNi    i§p£thi    TO    "'JELCl" 

^10  ^OR    i>l    TO    300 

50  EiiTER    i3Palhl  ;  V  i  (  I) ,  ZOt  I  ; ,  P3(  n  ,  vi  I  r,OiV  I  ;,'.)2\  r:,v3{  I  ) 

S8  NEXT    I 

70  U8=V^  L  > 

80  Ug  =  Li{  i  ; 

S0  FOR    I=Z    TO    300 

100  IF  i'i-;n>'jg  then  03=uicn 

118  IF  yisl)  03  THEN  '..'3  =  0in) 

123  NEXT  I 

130  O-{vS-U3)/I0 

US  FOR  1  =  1    "^0  300 

1B0  Z3=U8fD 

160  IF  ylc I  '  :Z3  THEN  X( 1  ;-0 

170  FOP  J=i  TO  9 

180  Zi=^S+; J-1 :-D 

190  Z2=^;3  +  J*D 

200  IF  ';i!l)\Zl  PiND  0]'I}--Z2  THEM  )((l)  =  d 

210  NEXT  J 

FOR  I-l  TO  300 

IF  K(  I  XS    THEN  SOTO  22S 

223  21=1/8  +  ;  10*0  ) 

224  Z2=0Sf; 3.7*0 i 

225  IF  ^l(I)'..Zl  oNO  01<n>22  THEN  X(I  )  =  0 
22B  NEXT  I 

230  GINIT 

240  !  PLOTTER  IS  705.'^HFGL" 

250  PLOTTER  IS  CRT ," INTERNAL" 

2G0  GRAPHICS  OM 

273  CSIZE  4.5, .£5 

280  MOUE  39,10 

2^0  Lf^BEL  "EMBEDDED  OORTEX  20  i1/5  " 

300  hOUE  4  3.5,7 

30 1  LABEL  -75%  FILM  COOLING" 

302  MOUE  41,85 

310  LABEL  "STREAMU'ISE  ^'ELOCITV 

320  MO'JE  56,17 

321  CSIZE  3.0 

330  LABEL  " Z   CM  " 

340  HOUE  5,5G 

341  CSIZE  3.0 

350  Le.be  1 5=  "V   Cr-1  " 

3S0  FOP  I-l  TO  B 

370      LABEL  LateUC  1,1] 
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390 
400 
•ilS 
4  20 
430 
4  40 


<-;-  s  '■■ 

4  4  S 
4  4? 

-t-  '1-  u 

■;  '  O 
<T  t  - 

4'3i 
4  -0  ^ 


455 
45c 


l/IEUPOR'i'    i5,ii4,:5,79 
SHQIJ    -ib,S.y,l;2 
FRnrlE 

Axes  .b,.5.0,3.;^z 

n  A  b  ';•     .  b  ,  .  'D  ,  -  1  D  ,  i  ^  .  _  ,  - 

CSI2E    :.3, .65 

iKX'E     -14.7,  10.55 

i. A B E L    ■'  './'El 0 C I T ■'■    R n r^^ 5 E  5 " 

A=10. 1 

FOR    K==i^    iU    9 

8~r-K " . 35 

M  f;  ( ,'  £    - 1  -I  .  3    9 

Z  l-==U3ff  I, -i-1  J  ♦D 

IF    K  =  9    THEN    Zj-'.'B^K*D-J-.7*D 

LABEL   'J3Ir.;G    '■D,2X,DC.D,X.DD.D- ;K,Z2,2i 

NEXT    K 

LABEL    USING    •■D,2X,D0.D,X,DD.D^' ;0,/l,V9 

C5IZE    1.3, .72 

F'OR    I-I    TO   300 

MOUE    Zea  /.YKI  ) 

IF    Yiax>S.0   AND    Z0':i)<-10.    iHE'^   GOTO   461 


4'JJ  3 

Ir      Yiv 1  'Jb. 

463 

LABEL  USING 

4Gi 

MEXT  I 

452 

CLIP  OFF 

4G3 

C5IZF  2.0, . 

FOR  1  =  -.:  TO  11.8  ;ST5F  2 


430 


LABEL  ■JSlilG    ■■:?  .OD.O-'i  II 


523  FOF  J--17.0  TO  5.0  STEP  2 

530      MOVE  J , - . ? 

548      Jl-J-H.O 

550       LABEL  US  IMG  "  Jt ,  r-^uO  .  0"  ;  J  1 

5b 0  NEXl  J 

570  £;>!D 
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18 

20 

30 
40 
50 
60 

.'  w 

f~»  *> 

b!5 

S0 

100 

110 

120 

133 

140 

158 

1S0 

170 

180 

130 

200 

210 

220 

230 

240 

250 

2G0 

270 

230 

290 

300 

310 

33S 
340 
350 
360 
370 
380 
330 
400 
410 
420 
430 
440 
450 
4B0 
470 
480 
4S0 
5S0 


RE.f1  FROGRnM  UECTOR 

REM  THIS  FROGRnil  U5E  Dr^Tn  !^RQM  *Jb[_C    DAIri  FILE! 

REfl  TO  PLOT  3LC0NDnRV  FLOW  VECTORS 

DIN  yi;e00  >.Z0(300 ',Py(800  ;,U<  30S  K'JUH^2  )  ,''2< 

Din  X^800),V{5),Z(5  ■ 

ASSIGN  ;aPath]  TO  "UEuCln'' 

FGR  1=1  TO  GO0 

ENTER  SFathl  ;  Via  ).ZC'-  I  ).F3':  I  ),U^  K'-'in  ).V2(  I 

NEXT  I 

GINIT 

Plotter  is  7?5,"hfsl" 

I  PLOTTER    IS    CRT,  ■■  liJTERNnL  " 

GRAPHICS  ON 

CSIZE  4.5, .GS 

HOUE  39, 3  0 

LABEL  " EMBEDDED  1'Ck''"F;<  l1  s '3  " 

LABEL    "75%    FiLri    COOL  IMG - 

riu'JE    35.  4.  85 

LABEL    '^SECONDARY    FLCw    UECTGRS" 

MOUE    G2, J  7 

CSIZE    3.0 

Lr.BEL     "Z       CM     " 

MOUE   5,G5 

CSIZE    5.0 

LabeIS="V      Cil    " 

FOR    1=1    ^0   G 

LABEL    LabeUI  1,1; 

r4EXT    I 
VIEWPORT    15,114,25,70 
SHOW    - IG , 5 , 0 , 12 
FRftnE 

r^XES    .5,  .5.0,0.2,2 
n  ,<  t  o     ,  -J  ,  .  b  ,  -  I  b  ,  1  ^  .  ^.  ,  ^ 
AXES    .5,. 5,5,12.2,2 
CSIZE    1.3, .72 
Sl=.2l 
A  1=0. 

FOR    K-7S1    TO    800 
Al=AU^l'vK  ) 
NEXT    K 
Al=Al*Si/40. 
MOUE    -14.5,10.3 
Y( 1  }  =  10.S 
Y(2)  =  10.B 
Z<  1  )  =  -14.5 
Z(2)=-14.5+nl*.l 
FOR    J=i    TO    2 

Plot  Z(  j  ),■/<] ; 

sEaT    j 

PFri'Ji^ 


>  3G0) 
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520 
S33 

54  Q 

550 
560 
570 

S80 
5S0 
G0G 
G10 


hO\iE    -15.  ,9.5 

C3IZE  2.0, .65 

LriBEL  "0.  1  F5  'JELOCITV" 

FOR  J-i  TO  BO0 

Yd  )  =  ri(  I  ) 

Z( 1  )-Z0(  I  ■ 

y(  2  )==v  1  i,  I  )  +  SI*U2\  I  ) 

Z( 2 --ZO? I  '-Sl-U3( I  ' 

I 

I  BUSINESS  TO  MAKE  r-PROWS 


■530 

Si'Zi 

G50 

GB0 

S70 

680 

S98 

700 

710 

720 

730 

740 

750 

7G0 

770 

780 

7G0 

800 

610 

820 

830 

840 

850 

8B0 

870 

880 

890 

900 

310 

920 

330 

340 

950 

9B0 

370 

S80 

990 

1000 

131.0 


Di-(  '.'(2  )-Y{  1  )  ) 
02-'  Z{  2  '-Z'  I  ;  ) 
UEG 
A3-r!TN{D2/Di  ; 

E3  =  3. 

ir  CI  :3  AND  D2.  Q 

iP  E3--1  .d    AND  A3 

04  =  ;  ^  Z;  2  >-2(  1  /  ;■■; 

[J4=.2G«04 

DS=i . 3 

D3-1.2 

IF  n4'-45  THEN  D3^ 


THEN  E3=l .0 
•■-45.  THEN  n4-A3+13t 

z+(  Y{2  /-Y(  1 )  yz )-  .5 


IF  Dl 
IF  Dl 
IF  Dl 
IF  Dl 
IF  Dl 
IF  01 
I'f 


IF  Dl 
IF  El 
IF  01 
IF  F2 

Y  (  3  )  = 
Z  V  3  )  = 

Y  (  4  )  = 
Z(4  i« 
A5=A4 
IF  01 
IF  Dl 
Y{5)  = 
2(5  )  = 


: 0  AND  02  ■  G 
>0  AND  02^3 
<0  AND  D2<0 
■0  AND  D2<:0 
<0  AND  02; 0 
<0  AND  A3>4 
■  --15.  AND  A 
>Q  AND  D2C0 
=1.0  AND  A3 
<0  AND  02 <0 
=1.0  AND  A3 


THEN  G9=-l .0 
THEN  D8=-l .0 
THEN  08=- i .0 


THEN 

D9  =  -i 

.0 

THEN 

D8  =  -l 

.0 

5  TriEi 

r^  09  =  1 

.9 

3-;0. 

THEN  0 

8^1 

.0 

THEN 

El=l. 

0 

.-45. 

THEN 

03  = 

1.0 

THEN  E2=i .0 


-45.  THEN  D9=l .0 
Y(2  )-D4*C0S( A4 )*D8 
Z(  2  )--D4*SINC  A4  )*D9 
Y(2> 
Z(2) 
-90 

.■0  AND  02  <0 
-:0  AND  D2<0 
Y<2  )+D4*CQS 
2(2  )  +  D4*SIN 


THEN  DS^l  .0 
THEN  09=1.0 

( A5  )>09 

(  AS  )*D8 


:nD  of  BUSINESS  TO  NAKE  ARROW; 


nOVE  20V  i  ;, Vi( I 
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1020  IF  Vl(I)-c.S  AND  Z0rn--10.  THLN  GOTO  1070 

1030  FOR  J=l  TO  G 

1040  PLOT  7(  J  ;,r(  J) 

1S50  NEXT  J 

iSGQ  PENUP 

1070  NEXT  I 

1080  CLIP  OFF 

1393  csize  :.0, .65 

1120  FOR  I--. 2  TC  II. £  STEP  2 

iliO  riG'.^b  -17.G,  I 

1120  J1-I+.2 

1130  LriBEL    USING     '?  ,  OD  ,  D"  ;  1 1 

1140  NEx'T    I 

1150  FOR  J--17.0  JQ   b.0  STEP  2 

1 150  ,^10t'E  J  ,  -  .  7 

1170  J].-=JiL.l5 

1180  LABEL  USING  " P , nDD. D" ; J  1 

liSe  NEXT  J 

1200  END 
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APPENDIX  C 
UNCERTAINTY  ANALYSIS 

Uncertainty  analysis  was  performed  using  the  method  originally  attributed  to 
Kline  and  McCormick,.  [Ref.  13].  Let  6^.  be  the  uncertainty  in  the  result  and  ^\  ^2  • 
.  .,  d^  be  the  uncertainties  associated  with  each  independent  variable.  The  uncertainty 
in  the  result  can  be  expressed  as 

\  =  [I((5R/^Xj)x(8.))]l''2  (eqnC.l) 

To  determine  the  uncertainty  of  the  pitch  angle,  a,  a  straight  line  approximation 
was  made  for  a 

a  =  m  +  b(Cpp-^^,j^)  (eqn  C.2) 

The  following  independent  variables  were  determined:  8^  =  ±1.97  o^  §^  =  ± 
2.32 o,  and  6^  =  ±  0.04637.  The  uncertainty  of  a,  was  determined  to  be  ±  2.36 ». 
The  high  uncertainy  of  the  pitch  angle  is  caused  by  the  probe  being  highly  sensitive  in 
the  pitch  plane. 

The  uncertainty  of  the  yaw  angle, P,  is  determined  in  a  similar  manner.  Once 
again,  a  straight  line  approximation  is  made. 

P  =  m  +  b(Cpy^^)  (eqn  C.3) 

In  this  case,  the  uncertainties  of  the  independent  variables  were  found  to  be  6^ 
=     ±    0.954,   6u    =    0.5924,   and  8(^„  =     ±    0.058.     From  these   values   the 

uncertainty  of  p  was  calculated  to  be  ±  1.29  <> 
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